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EXECUTIVE  PROGRAM  SUMMARY 


In  this  NRL  Memorandum  Report  we  describe  the  x-ray  laser  research  conducted 
by  the  Radiation  Hydrodynamics  Branch  under  BMDO  sponsorship  during  FY  1993. 
Some  of  the  experimental  work  was  performed  at  Sandia  National  Laboratories  in 
collaboration  with  NRL.  This  document  consists  of  three  sections  whose  authms  are 
separately  identified,  each  covering  in  detail  the  progress  in  a  different  aspect  of  the 
overall  effort 

The  work  of  FY  1993  has  brought  exciting  insights  and  discoveries  relevant  to  both 
principal  categories  of  x-ray  laser  drivers.  In  genial,  die  eneigizaa  of  x-ray  lasers  may 
be  divided  into  two  broad  groiq)s,  energy-rich  and  intensity-rich.  The  most  mergy-rich 
drivers  are  the  large  pulsed-power  genaatms,  which  can  dump  tens  of  kilojoules  or  more 
into  a  plasma  in  tens  of  nanoseconds,  via  an  electrical  discharge.  The  most  intensity-rich 
drivers  are  the  new  ultrashcat  pulse  lasers  which,  while  usually  coupling  a  joule  or  less 
to  the  medium,  do  it  in  a  fraction  of  a  picosecond,  over  a  small  focused  area  yielding 
intensities  of  1(F*  W  cm^*  or  higher.  Both  types  of  systems  have  specific  advantages.  A 
discharge  driver  is  able  to  produce  a  relatively  long  radiation  pulse  containing  considerable 
energy.  A  short-pulse  laser,  on  the  otho:  hand,  has  the  singular  capability  to  decouple 
the  ionization  state  of  the  plasma  fiom  its  temperature.  This  effect  can  occur  because 
field  ionization  dominates  the  usual  collisional  processes,  or,  alternatively,  because  the 
heating  of  the  electrons  outpaces  ionization  in  those  cases  where  collisional  effects  are 
important  In  either  instance,  a  uniquely  nonequilibrium  plasma  results  which  can  be 
highly  conducive  to  inversions  and  gain. 

For  several  years,  the  Na-Ne  resonantly  photopumped  x-ray  laser  has  been  pursued 
in  collaboration  with  Sandia  National  Laboratories  on  the  Saturn  facility  at  Sandia,  the 
world's  most  powerful  Z-pinch  driver.  This  effort  has  already  resulted  in  the  first  and  to 
date  the  only  demonstration  of  an  x-ray  population  invasion  driven  by  photopumping. 
During  FY  1993  a  noteworthy  event  in  this  continuing  research  was  the  first,  and  quite 

o 

strong,  detection  of  the  photopumped  4-2  He-like  line  complex  at  58  A.  No  attempt 
has  been  made  to  measure  gain  and  indeed,  this  wavelength  is  not  one  at  which  gain  is 
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expected  for  this  experiment.  However,  x-ray  laser  work  has  been  replete  with  instances 
of  both  unexpected  gain  and  unexpected  lack  of  gain,  and  in  many  cases  the  underlying 

o 

physical  causes  are  still  unresolved.  Relative  to  the  3-2  complex  at  78  A,  the  observed 
strength  of  this  line  is  at  least  a  factor  of  two  greater  than  that  predicted  by  either  the 
NRL  or  Lawrence  Livermore  National  Laboratory  models.  The  lines  which  are  expected 
to  show  die  greatest  gain  undor  the  Sandia  ^qierimental  conditions,  namely,  the  4-3 

o 

singlets  near  230  A,  have  not  yet  been  seen.  For  a  number  of  reasons  whose  details  are 
inappropriate  for  a  summary,  it  is  believed  that  these  lines  are  obscured  by  anomalously 
high  opacity  of  the  plastic  window  which  retains  the  neon  in  the  gas  cell.  Relatively  minor 
modifications  of  this  experiment  ought  to  be  able  to  more  effectively  use  the  radiation  of 
tile  Na  pinch  to  fully  ionize  this  window  and  tiiereby  reduce  its  abs(»ption. 

One  of  tile  most  attractive  methods  to  reduce  the  wavelength  and  increase  the 
efficiency  of  an  x-ray  laser  is  to  drive  an  inversion  relative  to  the  ground  state  of  an  atom 
or  ion.  For  instance,  for  the  Ly  a  line  of  a  hydrogenic  species,  the  transition  energy  is 
75%  of  the  ionization  potential,  compared  to  14%  for  the  3-2  line  which  can  lase  via 
recombination.  The  2-1  Ly  a  line  can  also  be  inverted  by  recombination  pumping,  but  it 
is  a  vdy  difficult  task  due  to  the  rapid  fillup  and  slow  emptying  of  the  ground  state.  This 
is  where  the  new  ultrashort  lasers  can  serve  a  critical  and  at  present  irreplaceable  role.  By 
stripping  an  ion  completely  and  rapidly,  these  devices  allow  the  development  of  transient 
inversions  which  could  lead  to  gain  in  Ly  a.  The  first  section  of  this  report  presents 
a  comprehensive  analysis  of  the  both  the  laser  specifications  and  plasma  conditions 
necessary  to  achieve  gain  in  the  Ly  a  lines  of  both  H  and  He.  It  is  found  that  for  these 
elements,  existing  systems  are  capable  of  demonstrating  such  a  laser. 

A  fundamental  analysis  of  the  interaction  of  such  short-pulse  lasers  with  a  target  is 
presented  in  the  second  section.  The  comprehensive  model  described  therein  includes 
detailed  atomic  physics  of  both  A1  and  Si,  non-LTE  radiation  hydrodynamics,  and  a 
Helmholtz  wave  solution  for  the  laser  field.  The  role  of  energetic  electrons  in  creating 
K  shell  vacancies  is  explored.  The  use  of  thin  multiple  layers  is  a  potentially  promising 
method  of  obtaining  information  about  ionization  and  (possibly)  also  the  electron  energy 
distribution.  In  some  of  the  upper  states  of  Li-like  Al,  inversions  are  calculated  to  develop 


v 


in  the  blowoff  plasma  for  the  higher  laser  intensities  considered. 

Tlie  final  section  demonstrates  dramatically  the  effect  that  decoupling  the  temperature 
and  ionization  state  in  the  presence  of  transient  population  behavior  can  have  on  the  gain  of 
a  well-known,  proven  laser  system.  In  the  case  of  Ne-Uke  Se,  gains  exceeding  600  cm“^ 
are  calculated  for  pump  laser  intensities  above  3  x  10^*  W  cm“^.  This  is  more  than  two 
ordears  of  magnitude  greater  than  present  expoimentally  observed  values.  Howev^,  given 
the  availability  of  superintense,  ultrashort-pulse  lasers,  there  are  sound  physical  reasons 
underlying  the  feasibility  of  such  a  huge  amplification.  Among  the  important  processes 
in  the  plasma  is  the  lagging  of  the  ionization  state  behind  the  electron  temp^ture.  In 
the  case  under  study,  Ne-like  Se  ions  are  exposed  to  coUisional  excitation  by  electrons 
much  hotter  than  would  occur  in  a  plasma  closo:  to  collisional-radiative  equilibrium.  This 
results  in  considerably  faster  excitation  of  the  iq)per  laser  levels,  and  is  responsible  in  part 
for  the  much  greater  predicted  gain. 
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Residual  energy  and  its  effect  on  gain  in  a  Lyman  a  laser 


P.  Pulsifer,  J  Jl  Apruzese,  J.  Davis  and  P.  Kepple 
Radiation  Hydrodynamics  Branch 
Hasma  Hiysics  Division 
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Abstract 


To  examine  prospects  for  gain  in  a  Lyman  a  recombination  laser  driven  by  a  high- 
intensity,  short-pulse  laser,  we  calculate  the  residual  energy  in  both  hydrogen  and  helium 
during  recombination  after  the  ionizing  pulse,  llie  expected  gain  as  a  function  of  residual  energy 
and  density  is  then  separately  evaluated.  The  residual  energy  calculation  includes  above-threshold 
ionization  (ATI)  in  the  presence  of  a  background  plasma,  as  well  as  inverse  bremsstrahlung 
heating.  At  electron  densities  over  10*^  cm~^  but  below  critical  density,  the  plasma  reduces  the 
ATI  energy  by  approximately  a  factor  of  two,  but  without  a  previously-reported  dependence  on 
the  pulsewidth.  Inverse-bremsstrahlung  heating  can  be  significant,  but  is  not  dominant  for  the 
parameters  considered.  Detailed  recombination-laser  gain  calculations  were  performed  for  the  Ly 
a  transitions  of  both  H  and  He,  using  Stark  profiles  to  represent  the  laser  line  cross  section.  To 
obtain  gain  of  near  2  cm~*  lasting  at  least  a  few  ps,  the  H  plasma  temperature  must  be  less  than 
3.5  eV  and  electron  density  between  4x  10*^  to  4x  10‘*  cm~^;  for  He,  the  temperature  must  be  less 
than  15  eV  and  electron  density  between  2x  10*^  and  2x  10'^  cm~^.  Our  calculations  indicate  that 
these  conditions  can  be  satisfied  for  H,  if  the  driving  laser  intensity  is  above  4  x  10'^  W  cm~^,  and 
for  He,  if  the  laser  intensity  is  above  1.7  x  10‘®  W  cm“^  and  wavelength  is  below  0.6  ftm. 
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1.  Introduction 

The  possibility  of  developing  a  recombination  x-ray  laser  which  is  driven  by  the  new  short- 
pulse,  high-intensity  table-top  lasers  has  recently  been  considered. With  sufficiently  high 
intensity,  a  population  inversion  can  be  obtained  with  large  transition  energies  (such  as  those 
involving  the  ground  state),  which  can  generate  significant  gain  at  UV  or  soft  x-ray  wavelengths. 
These  schemes  require  both  that  the  atoms  initially  be  ionized  above  the  lasing  stage  and  that 
the  residual  energy  of  the  electrons  be  low  to  facilitate  recombination  and  suppress  collisional 
excitation. 

In  this  paper,  we  analyze  the  conditions  required  for  gain  on  the  Lyman  a  transitions  of 
both  hydrogen  (1216  X)  and  helium  (304  X).  Such  lasers  were  initially  proposed  some  years 
ago.^  Our  gain  calculations  are  based  on  fiilly  self-consistent  Stark  profiles  for  the  completely 
stripped  background  gas.  The  residual  energy  is  calculated  including  the  electron  above-threshold- 
ionization  (ATI)  energy,  a  simple  model  of  plasma  effects,  and  inverse-bremsstrahlung,  or 
collisional,  heating. 

To  obtain  the  ATI  energy,  we  follow  a  well  established  procedure’*^  and  solve  for  the  classical 
(non-relativistic,  non-quantum)  electron  motion;  the  AH  energy  is  just  the  electron’s  energy  long 
after  the  pulse  has  passed  (after  several  ps,  when  recombination  is  occurring). 

PIC  simulations^  have  indicated  that  space-charge,  or  plasma  effects,  might  in  some  way 
significantly  enhance  or  diminish  the  eneigy  absorption  of  the  target  medium.  For  short  pulse 
lengths^  comparable  to  the  inverse  plasma  frequency,  space-charge  heating  apparently  stems 
from  the  non-adiabatic  ponderomotive  expulsion  of  electrons  (electron  cavitation),  which  induces 
plasma  oscillations.  An  earlier  study*  sought  to  model  plasma  effects  on  ATI  energy  by  including 
the  plasma  frequency  in  the  single-electron  equation  of  motion;  it  found  that  the  plasma-modified 
An  energy  was  fairly  large,  but  oscillated  by  several  orders  of  magnitude  as  a  function  of  electron 
density  and  pulse  period.  Here,  we  use  the  same  model,  but  our  findings  disagree  with  those 
results:  the  AH  energy  is  much  lower,  and  there  are  no  oscillations  with  electron  density.  For 
this  simple  model,  the  principal  effect  of  a  background  plasma  is  to  .substantially  increase  the 
residual  energy  near  the  critical  surface.  These  results  are  more  in  accord  with  those  of  other  recent 
theoretical  and  experimental  papers,^'**’  are  not  inconsistent  with  the  PIC  simulations  in  Ref.  8, 
and  are  in  fact  more  favorable  to  the  success  of  a  Ly  a  laser. 

One  potentially  important  factor  that  has  not  been  widely  considered  is  inverse  bremsstrahlung 
heating.^  Simple  scaling  arguments  indicate  that  this  could  be  a  much  larger  source  of  energy  to  the 
plasma  than  the  ATI  mechanism.  In  fact,  for  the  parameters  we  consider,  inverse  bremsstrahlung 
does  not  turn  out  to  be  the  most  important  source  of  heating,  largely  because  the  high  quiver  energy 
of  the  electrons  in  the  field  makes  the  electron-ion  collisional  cross  section  small  during  the  bulk 
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of  the  heating. 

The  organization  of  this  paper  is  as  follows.  In  Section  II,  we  consider  ionization  and  electron 
heating  by  the  laser  pulse,  and  show  how  to  obtain  the  residual  energy  during  recombination.  In 
Section  lU,  we  consider  the  recombination  process,  and  show  what  conditions  must  be  satisfied  by 
the  plasma  to  obtain  significant  gain.  These  results  will  show  the  optimal  pump  laser  parameters 
for  driving  the  recombination  laser.  Possible  experimental  parameters  are  discussed  in  the  final 
section. 


n.  Calculation  of  electron  residual  energy 


We  consider  a  laser  pulse  of  wavelength  A  (in  ftm)  and  peak  intensity  I  (in  W  cm~^),  incid  mt 
onto  an  initially  neutral  gas,  which  has  atomic  charge  Ze.  (An  electron  has  charge  — e  and  mass  m). 
The  laser  pulse  profile  is  assumed  to  be  such  that  the  electric  field  is  given  by  Bit)  =  E^nvit)  sinwt, 
with  the  pulse  envelope  given  by 


Eenv(t)=  Bo  scch 


2(t  - 
Tp 


) 


(1) 


and  the  peak  electric  field  Eq  =  y/lAQicI  (in  V/cm).  For  the  calculations  presented  here,  we 
consider  a  pulse  width  Tj,  =  100  fs  and  time  of  maximum  intensity  tma*  -  200  fs.  The  intense  laser 
field  multiply  ionizes  the  gas,  and  the  ion  densities  ni  (in  cm~^)  for  charge  state  Zi-i  evolve 
according  to  the  coupled  system  of  equations 

(2) 


where  Wi  is  the  rate  of  ionization  from  charge  state  Z,-i  to  Zi,  and  n.i  =  Wz-n  =  0.  The  evolving 
electron  density  is  Ug  =  Ziiii.  We  assume  that  the  laser  pulse  duration  is  much  shorter  than  the 
recombination  time,  so  that  during  the  pulse  recombination  can  be  neglected. 

There  are  various  theories  of  ionization  in  the  strong-field  regime.^  For  ATI  conditions, 
considered  here,  the  photon  energy  (fut)  =  1 .24/ A  eV)  is  less  than  the  ionization  energy  Ui,  which  in 
turn  is  less  than  or  equal  to  the  electron  quiver  energy  (defined  below).  Here,  tunneling  ionization 
might  be  expected  to  dominate,  and  the  ionization  rate  can  be  obtained  from  either  Keldysh"  or 
Ammosov  et  al.‘^  For  the  calculations  presented  here,  we  use  the  Ammosov  formula,*^  which  has 
some  experimental  support,’^  and  is  also  considerably  simpler  to  evaluate.  The  Ammosov  rate  is 


where  Wo.t,. 
atom: 


Wi(E)=  1.61  Wa.u. 


Z} 


n. 


*4.5 


10.87 


Z^  Ea 


T2n?-1.5 


nr  E 


exp 


2Z?  Eg 


3nr  E 


(3) 


=  4.1340X  10*®  sec“*,  nj  is  the  effective  principal  quantum  number  for  a  hydrogenic 


n,  = 


y/Ui/ 13.6  eV’ 


(4) 


A 


and  Ea.u.  -  5. 1421  x  10’  V/cm  is  the  electric  field  at  a  ground-state  hydrogen  electron. 

At  the  time  of  ionization,  we  assume  that  electrons  have  zero  energy;'  thereafter,  they  are 
(classically)  accelerated  by  fields  from  the  laser  and  from  the  surrounding  plasma.^  Most  of  the 
electron  energy  during  the  laser  pulse  comes  from  the  oscillatory  quiver  motion  in  the  laser  field, 
and  is  returned  to  the  field  when  the  pulse  passes.  The  relatively  small  residual  energy  after  the 
pulse  leaves  is  the  ATI  energy.  When  the  effect  of  the  plasma  is  ignored,  the  ATI  energy  is  the 
time-independent  part  of  the  energy,  which  is  determined  by  the  wave  phase  at  ionization.  >\fith  a 
plasma  oscillation  present,  there  is  no  time-independent  part  of  the  energy,  so  the  ATI  energy  must 
be  an  average  energy.  To  find  it,  we  follow  the  classical  trajectory  x{t)  of  an  electron,  with  the 
effect  of  the  charge  cloud  of  the  surrounding  plasma  being  treated  solely  as  a  spring-like  restraining 
force*  parameterized  by  the  plasma  frequency  Wp  = 

x+u)lx  =  -—E{t).  (5) 

^  m 

The  required  boundary  condition  on  Eq.  (5)  is  that  z  =  0  at  the  ionization  time  t  =  to.  >^th  the 
additional  imposed  condition  that  x(to)  =  0,  the  general  solution  to  Eq.  (5)  can  be  expressed  as 

[*  sinwpf  eE{t) _ ^ 


coswpf  /*  eE{i)  .  ^ 

z(<)  - —  I  - sin<«;pf at  ■ 

Wp  ^  m 


cosujpt  dt. 


If  the  electric  field  envelope  were  constant,  as  it  roughly  is  near  the  peak  of  the  pulse,  the  solution 
would  be  similar  to  that  obtained  elsewhere:* 

f  _  ^  r  ^  ^  Wpfo]  sin  u/pt 

fn(uj^  —  ufp)  I 

+  [— w  cosufto  sin Wpto  +  Up  sin uto  cos Upto\  cos Upt | .  (7) 

This  describes  the  complicated  electron  motion  during  the  pulse,  but  is  not  really  useful  for  finding 
the  ATI  energy,  which  is  detenr  ed  much  later.  The  electron  motion  at  those  late  times,  obtained 
from  the  <  — »  oo  limit  of  eq.  (6),  is  Just  simple  oscillation  at  the  plasma  frequency: 

r  *1 

= - [Asinu;pt  +  BcoswptJ  (8) 


A(to)=  f  sech  — — sinwfsinwpf 

Jto  L  J 

B(to)=  f  sech  — — sinutcosufptdt 

Jto  L  "^p  J 


2(f  -  tma.) 


sin  ut  cos  u)j,tdt 


The  time-dependence  of  the  factors  A  and  B  can  be  ignored  because  the  laser  pulse  (sech)  vanishes 
for  late  time;  they  depend  only  on  the  time  to  when  the  electrons  are  ionized.  If  the  plasma  period 
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is  smaller  than  the  times  of  interest  (the  ns  recombination  time),  the  single-electron  ATI  energy  is 
the  electron  energy  averaged  over  the  plasma-frequency  motion: 


c(fo)  = 


(II) 


where 

_  e^El 
’  4mu>^ 


(12) 


is  the  quiver,  or  ponderomotive,  energy. 

The  average  electron  residual  energy  (e)  after  the  passage  of  the  pulse  is  given  by  the  sum  of 
the  number  of  electrons  generated  at  each  to  times  the  residual  energy  of  these  electrons: 


E  J^ni(tQ)WjiE(to))e(to)dto 
_ 

T  S^niito)WiiE(to))dto 


(13) 


The  upper  limit  of  integration  is  chosen  as  Itnam  here  simply  for  convenience,  to  make  the 
integration  domain  symmetrical  over  the  pulse;  the  results  should  not  be  sensitive  to  the  value 
chosen  for  this  limit,  as  long  as  the  pulse  amplitude  is  small  there.  For  simplicity,  the  influence  of 
the  space-charge  field  on  the  ionization  rate  W(E)  has  been  neglected,  although  PIC  simulations^ 
have  indicated  that  electron  cavitation,  when  significant,  can  produce  a  space-charge  electric  field 
comparable  to  the  laser  field  E(t).  The  simple  model  used  here  could  be  extended  in  that  case 
by  deducing  the  space-charge  field  at  the  points  of  interest  (subject  to  details  of  pulse  profile  and 
ionization  dynamics)  and  then  using  the  total  field  in  eq.  (13).  We  define  the  residual  electron 
temperature  to  be  proportional  to  the  average  energy: 

r=|(e).  (14) 


When  the  laser  pulse  width  is  very  broad  compared  to  the  oscillator  frequency  and  plasma 
frequency  (i.e.,  wxp  >  1  and  WpXp  >1),  then  the  integrals  of  A  and  B  can  be  approximately 
evaluated,  and 


e(<o) 


e^E(to)^  Wp  sin^  u/fo  +  cos^  wfo 
4m  (u;^  —  u;|)^ 


(15) 


This  expression  shows  that,  long  after  the  pulse,  plasma  effects  do  not  decrease  the  residual  energy. 
There  is  no  periodic  dependence  on  the  plasma  frequency  in  the  residual  energy,  as  was  asserted 
elsewhere.*  That  periodic  variation  was  based  on  an  ATI  energy  calculated  from  some  average  of 
Eq.  (7),  and  so  incorrectly  included  electron  motion  during  the  pulse  in  the  ATI  energy  calculation. 
In  fact,  plasma  effects  can  increase  the  residual  energy,  but  the  effect  is  not  important  except  near 
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the  critical  surface,  where  the  plasma  and  wave  frequencies  are  near  resonance  and  the  residual 
energy  becomes  very  large. 

When  the  density  is  zero  (no  plasma)  the  single-particle  ATI  energy  from  Eq.  (15)  is 

c(f  o)  =  cos^  uitQ  ( 1 6) 


which  is  half  the  standard  zero-density  result.'  The  factor  of  two  difference  is  due  to  the 
averaging  over  a  plasma  period  performed  in  Eq.  (IS);  this  averaging  makes  sense  only  when 
the  plasma  frequency  is  large,  while  for  low  densities  the  peak  value  of  the  electron  energy 
should  be  used.  The  density  effect  on  residua!  energy  depends  on  the  ratio  =  ne/tio^, 
where  nern  1  12  x  lO^'/A^  is  the  critical  density,  where  the  plasma  becomes  opaque  to  the  laser 
radiation.  The  single-particle  ATI  energy  from  Eq.  (IS)  is  shown  in  Fig.  (1),  where  separate  curves 
are  drawn  for  different  ionization  phases  uiIq.  As  the  density  approaches  iia^,  the  ATI  energy 
sharply  increases,  and  becomes  less  dependent  on  the  time  of  ionization  Iq. 

Evidently,  plasma  does  not  significantly  increase  the  ATI  energy  at  much  less  than  O.lnertt. 
For  1.06  /tm  Nd:glass  laser  light,  the  ner«  =  10^’  cm~^,  while  for  a  248  nm  KrF  laser 
nerit  -2x  10^  cm~^;  thus,  according  to  our  simple  model,  plasma  modification  of  the  ATI  energy 
is  not  a  significant  factor  below  densities  of  10^  cm~^. 

Collisional  heating  from  inverse  bremsstrahlung  is  potentially  important,  possibly  dominant, 
in  dense  laser-heated  plasmas.^  The  inverse  bremsstrahlung  heating  rate  is  given  by 

=  y  «,  >'«  4  /o(0)  F(.vo/m)  (17) 

where  /o(v)  is  the  isotropic  part  of  the  electron  distribution  function,  vth  =  y/2kTlm  is  the  thermal 
velocity  (with  kT  the  electron  temperature  in  energy  units), 

i/.i  =  3.86  X  (18) 

is  the  electron-ion  collision  frequency,  with  Coulomb  logarithm  given  by 


log  A  =  22.4  +  log  [(Jfe7’)^/Vn*/2] , 


(19) 


and  F  is  a  correction  factor  to  the  cross  section  due  to  the  quiver  motion  of  the  electrons,'^ 
may  be  tq)proximated 


F(x)  = 


1 

1+ 2*2/3' 


which 

(20) 


The  rate  F/b  depends  on  the  electron  distribution  only  through  the  overall  normalization,  and  is  not 
sensitive  to  the  details  of  its  shape.  Therefore,  although  the  laser-produced  electrons  equilibrate 
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on  the  ps  to  tens  of  ps  timescale  (much  longer  than  the  pulse  width),  it  is  not  too  inaccurate  to 
approximate  the  distribution  function  to  be  a  Maxwellian,  so  that 

To  estimate  the  magnitude  of  inverse-'bremsstrahlung  heating,  we  determine  the  total  energy 
input  by  time-integrating  the  heating  rate,  Eq.  (17).  The  heating  rate  depends  on  the  time- varying 
electron  density,  obtained  by  solving  the  rate  equations,  Eq.  (2).  The  temperature  at  each  time 
is  given  by  the  average  energy  per  particle,  as  accumulated  from  ATI  (a  small  term),  prior 
inverse  bremsstrahlung  heating,  and  the  quiver  motion  (a  large  term).  The  inverse-bremsstrahlung 
contribution  to  the  final  temperature  is  just  the  total  average  energy  input  per  particle: 

Tib  =  z —  /  (Eib)  dt  (22) 

The  actual  residual  temperature  is  given  by  the  sum  of  contributions  from  ATI  (Eq.  (14))  and 
IB  heating  (Eq.  (22)).  Results  of  residual  energy  calculations  are  shown  in  Figs.  2-5.  The  contours 
in  Fig.  (2)  depict  the  average  degree  of  ionization  after  the  laser  pulse  of  an  initially  neutral  H 
gas,  as  a  function  of  laser  intensity  and  wavelength;  any  laser  power  above  the  Z  s  1  contour 
leaves  the  gas  fully  stripped.  The  contours  in  Fig.  (3)  show  the  residual  temperature  in  eV  of  the 
ionized  H  plasma.  The  inverse-bremsstrahlung  heating  rate  is  proportional  to  n.;  here,  an  atomic 
density  of  2x  10'*  cm~*  was  used,  where  there  is  relatively  high  gain.  Similarly,  Fig.  (4)  depicts 
the  average  degree  of  ionization  after  the  laser  pulse  of  an  initially  neutral  He  gas,  and  Fig.  (5) 
shows  the  residual  temperature  in  eV  as  a  function  of  laser  wavelength  and  intensity.  For  He,  the 
inverse-bremsstrahlung  heating  rate  was  computed  using  an  atomic  density  of  4  x  10'*  cm~’. 

m.  Conditions  required  for  gain  on  Ly  a  for  hydrogen  and  helium 

It  is  reasonable  to  assume  that  gain  on  Ly  a  via  transient  recombination^  would  be  easiest  to 
achieve  in  the  lowest  atomic  number  elements.  Therefore,  we  examine  H  and  He  closely  using 
time-dependent  multilevel  gain  calculations  to  establish  the  electron  temperature  (i.e.,  residual 
enei;gies)  and  densities  at  which  reasonable  amplification  would  occur.  These  can  be  related  to 
the  short-pulse  laser  parameters  which  are  predicted  to  provide  such  conditions  according  to  the 
calculations  of  Section  n.  In  this  way,  we  specify  what  would  constitute  a  promising  demonstration 
experiment. 

The  first  and  well-known  condition  for  transient  gain  to  the  ground  state  is  that  the  plasma  be 
initially  fully  stripped.^  Therefore,  laser  irradiances  must  be  used  which  are  at  or  above  those  in 
the  contour  plots  (Figs.  2-5)  corresponding  to  Z  =  1  for  H  and  Z  =  2  for  He.  In  all  the  calculations 
discussed  below,  the  plasma  is  assumed  initially  fiilly  stripped.  The  atomic  level  populations  are 
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calculated  as  a  function  of  time  for  a  fixed  assumed  density  and  electron  temperature.  The  ions  are 
assumed  to  remain  cold,  since  their  equilibration  time  (a  few  hundred  ps)  at  the  densities  of  interest 
far  exceeds  the  timescales  within  which  gain  is  expected.  The  algorithm  of  Ref.  IS  is  employed  to 
integrate  the  set  of  rate  equations,  viz., 

(23) 

i¥i  i¥i 

where  Ni  is  the  population  density  of  the  state  >  and  are  the  rates  connecting  all  states;  these 
atomic-level  rates  and  densities  are  not  to  be  confused  with  the  ionization-state  rates  and  densities 
of  Eq.  (2). 

In  Eq.  (23),  the  rates  connecting  the  various  levels  of  the  H  and  He  states  include  collisional 
excitation  and  de-excitation,  radiative,  three-body,  and  dielectronic  recombination,  collisional 
ionization,  and  radiative  decay.  Most  of  the  atomic  data  for  these  simple  and  well-studied  species 
are  readily  available  from  numerous  sources  in  the  literature.  Other  rates  were  calculated  using 
standard  techniques  (e.g.,  the  semiclassical  approximation  for  collisional  rates  corresponding  to 
radiatively  allowed  collisions).  The  levels  of  the  one-electron  species  are  distinguished  only  by 
principal  quantum  number  n,  although  fine  structure  effects  on  Ly  a  are  included  in  the  overall 
Stark  profile  used  in  calculating  gain. 

The  highest  principal  quantum  number  included  in  the  atomic  model  is  n  s  5  for  both  H  and 
He.  This  is  tqiproximately  in  accord  with  the  loglis-Teller  criterion  as  updated  by  Griem  (Ref.  16, 
p.  12S):  At  an  electron  density  of  10‘*  cm~^.  Ref.  16  predicts  that  the  series  limit  occurs  at  n  s  4-S 
for  H,  and  n  s  6-7  for  He.  While  we  do  not  include  a  self-consistent  treatment  of  continuum 
lowering,  it  is  clearly  preferable  to  truncate  the  levels  at  n  =  5  rather  than  at  some  arbitrary  higher 
number  (10  or  20,  for  instance),  far  above  the  last  discrete  level  at  densities  of  interest. 

The  gain  g  in  cm~*  on  the  n  =  2  to  1  transition  of  a  hydrogenic  species  is  given  by 

g  =  N2<r2i  -  Niai2  (24) 

where  0-21  (a^)  are  the  emission  (absorption)  cross  sections  in  cm^  for  the  transition  at  the  peak 
of  the  profile,  and  the  N's  are  the  population  densities  of  the  n  =  1  and  n  =  2  levels.  In  Eq.  (24) 
a  positive  number  corresponds  to  amplification  and  a  negative  one  to  absorption.  For  Ly  a, 
<r2i  =  <7'i2/4,  and  n  =  2  must  have  at  least  4  times  the  population  of  the  ground  state  to  obtain  gain. 
Clearly,  this  can  only  occur  on  a  transient  basis  since  population  tends  to  pile  up  in  n  =  1,  which 
has  no  radiative  decay  channel.  However,  three-body  recombination  from  the  bare  nucleus,  per 
unit  statistical  weight,  is  approximately  proportional  to  n*/T.  Therefore,  one  anticipates  a  brief 
period  of  gain  for  an  {q)propriate  range  of  sufficiently  low  temperatures,  if  n  =  1  is  initially  empty. 

A  key  feature  of  the  present  work  is  that  we  have  employed  the  best  available  representation 
of  the  Ly  a  cross  section,  namely  self-consistently  calculated  Stark  profiles  for  both  H  and  He.’^ 
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For  H,  we  find  that  a-^  is  approximately  given  by 

/|0l8\  0  *9 

(25a) 

whereas,  for  He: 

=  (256, 

For  T  <x  Z^,  the  He  cross  section  exceeds  fiiat  of  H  by  a  factor  of  4  at  tie  =  10‘*  cm“^. 
Partly  offsetting  this,  however,  is  the  fact  that  there  are  twice  as  many  H  ions  per  electron  as  He 
ions.  Therefore,  we  expect  gains  of  comparable  magnitude.  Also,  note  from  Eq.  (25)  that  as  the 
density  increases,  the  Staric  broadening  results  in  a  nearly  inversely  proportional  decrease  of  the 
peak  cross  section.  However,  gain  can  still  increase  with  density  due  to  the  increase  in  three-body 
recombination  and  the  greater  number  of  ions  available  for  stimulated  emission. 

A  subset  of  the  calculations  which  have  been  carried  out  for  H  and  He  is  presented  in 
Figs.  6-10.  There  is  no  lower  temperature  limit  below  which  gain  is  not  possible;  indeed,  due  to 
the  T~*  dependence  of  three-body  recombination,  gain  increases  at  lower  temperatures,  assuming 
of  course  that  (he  plasma  is  initially  fully  stripped.  To  obtain  a  gain  of  ~  2  cm~‘  lasti;  g  at  least  a 
few  ps,  the  electron  temperature  in  a  He  plasma  must  be  <  15  eV,  whereas  for  H  the  electrons  must 
be  no  hotter  than  3  J  eV.  An  upper  limit  on  the  density  is  set  by  the  desired  gain  duration.  Note 
from  Figs.  6-10  that,  as  expected,  the  gain  increases  wifii  density,  but  decreases  with  increasing 
temperate.  However,  the  duration  of  gain  decreases  sharply  with  increased  density  since  the 
absorbing  n  s  1  level  fills  up  correspondingly  more  rapidly.  For  gain  duration  of  at  least  a  few  ps, 
the  electron  density  must  be  no  higher  than  ~  2  x  10*’  cm“^  for  He,  or  4  x  10**  cm~*  for  H.  Since 
the  lifetime  of  the  H  Ly  a  transition  (n-averaged)  is  2.1  ns,  compared  to  133  ps  for  He,  one  might 
initially  expect  considerably  greater  gain  duration  for  H.  However,  the  corresponding  collisional 
mixing  rates  are  nearly  an  order  of  magnitude  faster  for  H  at  a  given  density,  and  the  interplay  of 
these  effects  results  in  comparable  gain  durations  for  both  elements  (ranging  from  a  few  ps  to  tens 
of  ps). 

Refraction  of  the  x-ray  laser  beam  can  pose  a  serious  threat  to  the  practical  attainment  of 
amplification,  due  to  the  density  gradients  expected  in  short-pulse  laser-driven  plasmas.  This  issue 
has  been  analyzed  by  Amendt  et  al.,^  who  assumed  a  Gaussian  beam  profile.  If  6  is  the  beam  waist 
radius,  the  x-ray  lasing  length  Zjt  is  approximately  given  by 

Zji  =  1.26^^  (26) 

where  na^  is  the  critical  density,  defined  above.  Typical  electron  densities  for  good  gain  in  H  and 
He  are  2  and  4  x  10**  cm“*,  respectively  (cf.  Figs.  6-10).  We  take  the  beam  waist  radius  to  be 
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13  fim,  as  measured  in  a  recent  demonstration  of  optical  guiding  of  intense  laser  beam  pulses.'^ 
For  these  conditions,  Eq.  (26)  gives  a  propagation  length  of  0.3  cm  for  H  and  0.9  cm  for  the 
considerably  shorter  wavelength  of  the  He  laser.  Figs.  8  and  9  indicate  that  gain  length  products 
(gl)  greater  than  or  equal  to  4  are  achievable  for  He  if  T*  <  5  eV.  For  H,  Te  <  1  eV  is  necessary  to 
achieve  a  similar  gl.  As  pointed  out  in  Ref.  5,  a  larger  beam  waist  radius  would  significantly  relax 
these  requirements.  In  addition,  the  presence  of  a  beam  guiding  channel’^  with  a  more  favorable 
refractive  index  profile  would  very  likely  lead  to  considerably  longer  propagation  lengths. 

rv.  Summary  and  conclusions 

We  have  presented  a  procedure  to  calculate  the  residual  energy  of  a  plasma  during  the 
recombination  period  after  passage  of  a  high-intensity,  short  laser  pulse.  We  then  calculated  the 
expected  residual  energy  and  degree  of  ionization  for  H  and  He  over  a  range  of  wavelengths  and 
laser  intensities.  These  results  are  shown  as  contour  plots  in  Figs.  2-5.  The  inclusion  of  a  plasma 
oscillation  does  not  significantly  modify  the  residual  electron  energy  (other  than  the  factor  of  two 
reduction  due  to  the  residual  plasma  oscillation,  as  discussed  at  eq.  (16)),  as  long  as  the  plasma 
density  is  not  close  to  the  critical  density. 

While  our  model  does  not  include  it,  we  do  not  expect  the  space-charge  heating  reported  in 
PIC  simulations^*'*’  to  be  important  for  parameters  where  significant  gain  should  be  found.  The 
short  pulse  width  limits  heating  due  to  instability  growth  and  variations  in  drift  velocities.'*’  The 
laser  power  required  is  not  ultra-high  (i.e.,  I  <  10'*  W  cm“^),  and  the  density  is  relatively  high,  so 
both  electron  cavitation  and  its  effect  on  the  final  plasma  temperature*  should  be  minimal;  even  if 
this  were  not  true,  electron-cavitation  heating  seems  to  be  much  less  important  when  ufpTpflic  >  1 , 
which  is  the  case  for  the  high-gain  densities  given  here. 

Although  certainly  not  conclusive,  it  is  gr^fying  to  note  that  our  results  are  close  to  recent 
experimental  measurements  by  Mohideen  et  al.'^  Their  experiment  measured  residual  energy 
in  a  low-density  (7  x  10*  atoms  cm"’)  He  gas,  where  the  factor  of  1/2  from  averaging  over  the 
late-time  plasma  oscillation  in  eq.  (16)  should  be  omitted.  The  laser  wavelength  was  0.820  /im, 
with  intensity  7  x  10'’  W  cm"^  and  pulsewidth  180  fs.  Although  their  pulsewidth  was  slightly 
longer  than  assumed  here  (which  would  tend  to  increase  the  residual  energy).  Fig.  (S)  shows  that 
for  these  parameters  the  predicted  residual  temperature  is  roughly  20  eV,  compared  to  the  30  eV 
measured  best-fit  temperature. 

Through  the  residual  energy  and  ionization  calculations,  we  have  established  the  plasma 
conditions  during  recombination.  We  have  then  used  an  inclusive  atomic  model  to  calculate  gain 
and  those  conditions  for  which  significant  gain  is  obtained.  These  results  are  shown  in  Figs.  6-10. 
The  region  in  the  contour  plots  where  the  plasma  is  both  fully  stripped  and  has  residual  energy 
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within  the  given  limits  is  where  gain  is  expected.  For  H,  fully  stripped  plasma  is  obtained  at 
intensities  above  4 x  10*'*  W  cm~^;  the  necessary  electron  temperature,  less  than  3.5  eV,  is  satisfied 
at  all  wavelengths.  For  He,  fully  stripped  plasma  is  obtained  when  the  laser  intensity  is  above 
1.7  X  10‘*  W  cm~^;  if  this  is  true,  the  electron  temperature  is  below  the  maximum  15  eV  when  the 
wavelength  is  below  0.6  ^m. 

The  minimum  density  for  which  adequate  gain  is  obtained  can  be  seen  from  Figs.  6-10  to  be 
roughly  10*’  cm~^.  Although  the  residual  energy  increases  with  increased  density,  mostly  because 
of  inverse  bremsstrahlung,  the  most  important  effect  is  the  decrease  in  gain  duration  (for  a  fixed 
residual  energy)  with  increased  density;  this  limits  the  electron  density  to  roughly  4x  10*’  cm~^ 
for  H,  or  2  X 10*’  cm“’  for  He. 

The  simple  model  used  here  notably  omits  consideration  of  collisions.  At  typical  laser-heated 
plasma  parameters  of  ibT=10  eV  and  n  s  2x  10*’  cm~’,  the  collision  time  is  about  0.8  ps.  Thus, 
these  plasmas  are  marginally  coUisional,  and  collision  effects  should  be  studied.  Still,  collisions 
do  not  clearly  dominate,  and  there  is  probably  some  validity,  at  least  for  simple  estimates,  to  the 
residual  electron  energies  given  here. 
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Figure  Captions 

Figure  1.  Effect  of  the  plasma  oscillation  on  the  single-particle  ATI  energy.  The  ATI  energy,  in 
units  of  the  quiver  energy  e,,  is  shown  as  a  function  of  density,  in  units  of  the  critical  density 
Aerit,  for  different  times  of  ionization  within  the  pulse  wto-  The  solid  line  is  for  wto  =  0,  when  the 
zero-density  ATI  energy  is  maximum;  the  dashed  line  is  for  uto  =  x-/4;  and  the  dotted  line  is  for 
<irfo  =  ir/2,  when  the  zero-density  ATI  energy  is  zero. 

Hgure  2.  Residual  ionization  contours  for  initially  neutral  H.  Contours  represent  average  ionization 
level  Z  -  nefriH'f  for  laser  intensity  above  the  Z  =  1  contour,  the  gas  is  left  fully  stripped. 

Figure  3.  Contours  of  residual  temperature  in  eV  of  H,  as  described  in  Fig.  (2).  For  inverse- 
bremsstrahlung  heating  calculations,  an  initial  atomic  density  of  2  x  10*^  cm~^  is  assumed. 

Figure  4.  Residual  ionization  contours  for  initially  neutral  He.  Contours  represent  average 
ionization  level  Z  =  n^InHe,  for  laser  intensity  above  the  Z  =  2  contour,  the  gas  is  left  fully 
stripped. 

Figure  S.  Contours  of  residual  temperature  in  eV  of  He,  as  described  in  Fig.  (4).  For  inverse- 
bremsstrahlung  heating  calculations,  an  initial  atomic  density  of  4  x  10'^  cm~^  is  assumed. 

Figure  6.  Ly  a  gain  coefficient  for  H  as  a  function  of  time,  assuming  an  initially  fully-stripped 
plasma  of  the  indicated  atomic  densities.  Electron  temperature  is  1  eV;  ions  are  assumed  cold,  and 
Stark  profiles  are  employed  to  obtain  the  gain  cross  section. 

Figure  7.  As  in  Fig.  6  except  that  the  electron  temperature  is  2  eV. 

Figure  8.  Ly  a  gain  coefficient  for  He  as  a  function  of  time,  assuming  an  initially  fully-stripped 
plasma  of  the  indicated  atomic  densities.  Electron  temperature  is  2  eV;  ions  are  assumed  cold,  and 
Stark  profiles  are  employed  to  obtain  the  gain  cross  section. 

Figure  9.  As  in  Fig.  8  except  that  the  electron  temperature  is  5  eV. 

Figure  10.  As  in  Figs.  8  and  9,  except  that  the  elytron  temperature  is  10  eV. 
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Ultra  Short  Laser  Produced  Al/Si  Plasma 


J.  Davis,  R.  Clark  and  J.  Giuliani 
Plasma  Physics  Division 


ABSTRACT 

The  x-ray  emission  from  an  ultra-short  pulse  laser  produced  aluminum/silicon  plasma  is  investigated. 
A  time  history  of  the  L-  and  K-shell  line  spectra,  including  the  Ka  lines  and  continuum  radiation  is 
presented  for  a  laser  intensity  of  1(P^  watts/cm^.  The  simulations  are  performed  using  a  non-LTE  radiation- 
hydrodynamics  model  coupled  to  a  Helmholtz  wave  equation  for  P-wave  polarization.  A  fraction  of  the 
abstnbed  oiergy  is  expended  in  the  production  of  fast  electrons,  which  are  transported  and  deposited  into 
the  cold  material.  These  electrons  create  K-shell  vacancies  which  produce  characteristic  Ka  line  radiation. 
Analysis  of  this  radiation  provides  information  about  the  ionization  of  the  target  and  the  spectrum  of  the 
energetic  electrons. 

1.  INTRODUCTION 

The  interaction  of  high  powered  pulsed  lasers  with  planar  targets  has  been  the  subject  of  intensive 
investigation  for  over  three  decades.  The  major  enq>hasis  of  the  research  was  focused  on  stabilizing  the 
Rayleigh-Taylor  (RT)  instability  for  application  to  laser  fusion  experiments.  Although  interest  remains 
hig6  for  the  eventual  control  of  the  RT  instability,  a  necessary  criterion  for  the  success  of  the  laser  fusion 
program,  a  number  of  other  activities  and  technologies  has  emerged.  One  of  these  activities  results  from 
the  development  of  techniques  which  makes  it  possible  to  produce  ultra  short  pulse  lasers  with  pulsewidths 
of  tens  to  hundreds  of  femtoseconds  and  intensities  in  excess  of  10^'  watts/cm^.  This  newly  evolving 
technology  is  creating  a  number  of  innovative  and  exciting  applications  involving  the  interaction  of  ultra 
short  pulse  lasns  with  matter  including  x-ray  lasers,  laser  produced  solid  density  plasmas,  subnanosecond 
spectroscopy,  and  the  interaction  of  intense  laser  fields  with  atoms,  molecules  and  energetic  electron  beams. 
The  intraaction  of  these  intense  lasers  with  planar  targets  makes  it  feasible  to  investigate  solid  state  plasmas 
as  a  potential  source  of  intense  x-ray  radiation  in  a  relatively  quiescent  environment,  since  the  target  plasma 
remains  nearly  stationary  during  the  initial  burst  of  x-rays. 

In  this  paper,  we  investigate  the  interaction  of  an  ultra  short  pulse  laser  with  a  planar  aluminum  target, 
with  the  focus  of  the  research  directed  toward  characterizing  the  radiative  properties  and  behavior  of  the 
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evolving  plasma.  The  resulting  radiation  pulse  represents  a  potentially  useful  source  of  x-rays  for  a  variety 
of  applications.  In  addition  to  the  thermal  plasma  reservoir,  there  are  groups  of  energetic  electrons  created 
in  the  absorption  process  that  deposit  their  energy  in  the  colder  dense  regions  of  the  target  producing  inner 
shell  transitions  on  a  time  scale  of  a  few  picoseconds.  The  mechanisms  for  the  production  and  deposition 
of  energetic  electrons  have  recently  been  studied  (Chen,  Soom,  Yaakobi,  Uchida  and  Meyerhofer  1993)  by 
making  use  of  the  resulting  iimer  shell  radiation.  For  example,  a  series  of  Ko  lines  is  produced  by  resonance 
fluorescence  in  various  ions  of  the  target  material.  The  relative  strengths  of  these  lines  can  be  used  to 
estimate  ionic  abundances  in  the  vicinity  of  the  energetic  electron  deposition  as  well  as  the  characteristics 
of  the  fast  electrons. 


2.  MODEL 

2.1  Hydrodynamic  Model 

The  basic  hydrodynamic  variables  consisting  of  mass,  momentum  and  total  energy  are  transported  in 
one  dimension  using  a  Lagrangian  flux-corrected  transport  numerical  scheme.  The  hydrodynamic  equations 
solved  are 
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where  p  is  mass  density,  u  is  velocity,  P  is  pressure  and  ep  is  total  energy,  irad,  iam,  cia«er  and  are 
the  rates  of  energy  gain  (or  loss)  due  to  radiation,  thermal  conduction,  laser  deposition  and  fast  electron 
deposition,  respectively.  The  ponderomotive  force,  Fp,  is  given  by 
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where  u  is  the  incident  laser  frequency,  upe  is  the  electron  plasma  frequency  and  E  is  the  local  electric 
field. 

A  single  temperature  approximation  is  employed,  T*  =  T<,  where 
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where  n,  is  the  electron  density,  e,  is  the  ionization  energy,  is  the  ion  density  of  species  m,  and  Pp  is  the 
Fermi  pressure.  The  assumption  of  a  single  temperature  breaks  down  in  situations  where  the  equilibration 
time  is  comparable  with  the  heating  time  and  the  local  average  charge  state  is  small.  Near  the  critical 
surface,  for  example,  the  equilibration  time  condition  is  violated  long  before  the  peak  of  the  laser  pulse, 
but  because  of  the  local  ionization  (initially  from  multiphoton  ionization),  an  acceptably  small  error  in  the 
electron  tempenUure  is  incurred. 


22  Atomic  Model 


The  ionic  populations  in  the  plasma  may  be  characterized  by  a  set  of  atomic  rate  equations  of  the  form 

-  Wuk„.fuu„), 


where  /am  is  the  fractional  population  of  atomic  level  i  for  species  m  in  spatial  zone  k,  and  TT^am  is  the  net 
rate  describing  the  transition  from  initial  state  j  to  final  state  i.  An  equation  of  this  type  is  constructed  for  each 
of  the  atomic  levels  included  in  the  model.  The  rate  coefficients  that  are  used  to  calculate  the  populating 
and  depopulating  rates,  are  calculated  using  various  scattering  techniques  and  the  methods  used 

in  calculating  the  corresponding  rate  coefficients  are  summarized  in  Duston,  Qark,  Davis  and  Apruzese 
(1983). 

Radiation  emission  from  the  plasma  and  its  opacity  are  dependent  on  the  local  atomic-level  population 
densities.  Except  for  optically  thin  plasmas,  however,  the  level  populations  depend  on  the  radiation 
field,  since  optical  pumping  via  photoionization  and  photoexcitation  can  produce  significant  population 
redistribution.  Thus,  the  ionization  and  radiation  transport  processes  are  strongly  coupled  and  must  be 
solved  self-consistently.  In  this  model,  an  iterative  procedure  (Duston  et  al.  1983)  is  used,  where  level 
populations  are  calculated  using  the  radiation  field  from  the  previous  iteration,  then  using  these  populations 
to  calculate  a  new  radiation  field  until  convergence  is  reached. 

The  prescription  for  multiphoton  ionization  follows  the  tunneling  ionization  theory  of  Anunosov, 
Delone  and  Krainov  (1986),  given  by 
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where  =  4.1x10“*"^  is  the  atomic  unit  of  frequency,  B,.,.  =  5.lxlO®V/cm  is  the  atomic  field  strength, 
and  the  effective  quantum  number  ng//  is  obtained  by  hydrogenic  scaling 
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where  Uh  =  iZ.fieV  is  the  ionization  potential  of  hydrogen,  Ui  is  the  ionization  potential  of  the  critical 
electron  orbit,  and  Z  is  the  residual  charge  seen  by  the  electron.  Multiphoton  processes  produce  some 
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initial  ionization  of  the  plasma,  but  become  unimportant  compared  with  coUisional  ionization  as  the  plasma 
is  heated. 

The  atomic  model  for  aluminum  used  in  this  calculation  consists  of  the  ground  states  and  1 18  excited 
levels  through  neutral  aluminum,  including  structure  up  to  N=:10  for  H-like,  and  to  N=9  for  He-like  and  U- 
like  aluminum.  The  atomic  model  for  silicon  consists  of  the  ground  states  and  65  excited  levels.  Ionization 
lowering  is  accounted  for  by  means  of  an  ion  sphere  model;  the  bound-free  radiation  is  limited  from  states 
which  merge  with  the  continuum.  Radiation  transport  is  carried  out  using  a  probability-of-escape  formalism 
which  is  described  in  Apruzese  (1981  and  1985). 


23  Laser  deposition  model 


Hie  absorption  of  a  laser  pulse  obliquely  incident  on  a  planar  target  is  calculated  by  solving  the 
Helmholtz  wave  equation.  In  this  way,  the  incident,  reflected  and  transmitted  waves  are  self-consistently 
characterized  in  each  spatial  zone  in  the  medium.  The  conductivity  is  obtained  from  the  Cauble  et.  al. 
(1987  and  1990)  model.  The  conductivity  is  complex,  with 

» =  — 
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where  p  is  the  collision  frequency  and  w  is  the  laser  frequency. 

The  Helmholtz  equations  for  the  electric  and  magnetic  fields  can  be  written 

V*B  +  fl  +  i  — )  E  =  V’B  +  ^  e'  B  =  0 
c*  \  w  /  c* 

V’B  -  ^x(VxB)  +  ic'B  =  0 

where  e'  is  the  plasma  dielectric  function.  For  S-polarization,  take  E  =  B,(*,y)e,  =  In 

this  case,  it  is  convenient  to  use  the  Helmholtz  equation  for  E,  which  reduces  to 
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For  P-wave  polarization,  it  is  mote  convenient  to  solve  for  the  magnetic  field,  B  =  and  the 

corresponding  Helmholtz  equation  can  be  written 
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A  Runga-Kutta  numerical  scheme  is  employed  to  solve  the  appropriate  Helmholtz  equation.  An 
adiq)tive  gridding  is  used  to  resolve  the  fields  in  the  vicinity  of  the  critical  surface.  The  computational 
mesh  used  for  the  solution  of  the  fields  is  much  finer  than  that  used  for  the  hydrodynamics  and  radiation 
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transport;  logarithmic  interpolation  is  employed  to  obtain  the  electron  density  and  collision  frequency  from 
the  hydrodynamic  computational  mesh. 

2.4  Energetic  electron  transport 


A  fixed  fraction  of  the  Uier  energy  deposited  in  each  spatial  zone  (10  percent  in  the  present  calculation) 
is  assumed  to  be  expended  in  the  production  of  energetic  electrons.  These  fast  electrons  are  assumed  to 
be  monoenergetic  (at  an  energy  between  4.0  and  10.0  keV).  This  assumption  is  made  to  simplify  the  the 
interx»etation  of  the  resulting  data;  a  Maxwellian  (or  an  arbitrary  energy  distribution)  could  be  enq>loyed 
in  the  model.  The  energetic  electrons  are  transported,  and  subsequently  deposit  their  energy  in  the  plasma, 
giving  rise  to  local  heating  and  irmer-shell  vacancies. 

At  sufficiently  high  electron  energies,  the  Bethe  theory  is  reasonably  accurate.  Then  the  fast  electron 
stopping  power  can  be  written  in  the  Bethe-Bloch  form  (as  given  in  Powell  1976) 


8x 


2Te*ni 

B 


[Zb  ln(^^)] 


i.ieEx 


where  x<n  is  an  average  ionization  potential,  Xo  is  an  effective  potential  related  to  the  electron  plasma 
frequency,  and  Zb  and  are  the  effective  number  of  bound  and  free  electrons  per  ion.  The  continuous 
slowing  down  model  breaks  down  when  quantum  transitions  are  of  the  same  order  in  energy  as  the  beant 


2.5  Ka  radiatioii  model 

The  K-shell  vacancies  which  are  created  are  r^idly  filled  by  either  Auger  decay  (radiationless  cascade) 
or  resonance  fluorescence  (valence  electron  fills  vacancy  with  the  emission  of  Ka  radiation).  The  energy- 
dependent  Ko  production  rate  is  calculated  for  each  aluminum  ionization  stage  in  each  spatial  zone,  from 
A1 1  (neutral  aluminum)  through  A1  DC  (boron-like  aluminum),  and  for  the  corresponding  silicon  ionization 
stages.  The  K-shell  vacancy  cross  sections  for  aluminum  and  silicon  were  obtained  from  the  calculations  of 
Blaha  (1992),  and  represent  a  small  fraction  of  the  total  stopping  power.  From  threshold  to  about  lOkeV,  the 
energy  dependence  of  these  cross  sections  can  be  rq>proximated  by  a  simple  slowing-down  approximation 

a(B)  =  7^--— r  1“(— ) 

(E  Xk)  'xk' 

where  C  for  aluminum  is  equal  to  9.5xl0~^*cm’/eV'’  and  xk  is  given  in  Table  I  for  each  ionized  state  of 
aluminum.  The  fluorescence/Auger  decay  branching  ratios  (taken  from  Duston  et  al.  1983)  vary  with 
ionization  stage,  and  are  typically  about  5  percent.  The  Ka  line  widths  are  determined  from  the  inverse 
lifetimes.  The  Ka  lines  for  Al  I-IV  (and  for  Si  I-V)  lie  so  close  together  in  energy  that  they  appear 
qiecttoscopically  as  a  single  line,  whereas  those  of  Al  V-DC  (and  Si  VI-X)  are  sufficiently  separated  that 
th^  can  be  resolved. 
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Table  I  lists  by  ion  the  Ka  energy,  the  resonance  fiuorescence/Auger  decay  branching  ratio,  and  the 
potential  used  to  fit  the  K-shell  vacancy  cross  section  for  aluminum.  A  similar  set  of  energies  and  branching 
ratios  is  used  for  the  silicon  radiation. 

TABLE  I:  Kq  Parameters  for  Aluminum 


Ion 

Ka  Energy 

Ka /vacancy 

Xft 

A1 

I 

1486.99  eV 

0.0330 

1600  eV 

A1 

II 

1487.52  eV 

0.0460 

1630  cV 

A1 

m 

1488.24  eV 

0.0450 

1660  eV 

A1 

IV 

1488.95  eV 

0.0434 

1690  eV 

A1 

V 

1499.58  eV 

0.0479 

1750  eV 

A1 

VI 

1511.64  eV 

0.0522 

1810  eV 

A1 

vn 

1525.40  eV 

0.0619 

1870  eV 

A1 

VIII 

1540.37  eV 

0.0700 

1930  cV 

A1 

IX 

1557.01  eV 

0.0623 

2000  eV 

3.  RESULTS 

3.1  Alttminum  Ihrget 

An  ultrashort-pulse  laser  (0.22  psec  FWHM  duration  Gaussian  with  no  prepulse)  was  assumed  to  be 
obliquely  incident  (30  degrees  from  the  normal)  on  solid  density  aluminum  and  layered  aluminum/silicon 
targets.  The  peak  power  was  1.0x10^^  W/cm^  at  a  wavelength  of  0.26  Figure  1  shows  the  time  dependence 
of  the  incident  laser  power,  the  absorption  fraction  (absorbed  energy  relative  to  the  total  incident  energy) 
a.ssuming  P-wave  deposition,  and  %e  temperature  at  the  critical  surface. 

The  reflectivity  varies  from  about  40  to  70  percent  during  the  laser  pulse,  and  is  about  50  percent  near 
peak  incident  power.  In  the  case  of  layered  aluminum/silicon  targets,  the  details  of  the  absorption  may  be 
different,  but  the  laser  deposition  features  which  are  important  to  this  discussion  (reflectivity,  absorption 
fraction,  etc.)  will  be  approximately  the  same.  The  energy  distribution  history  is  plotted  in  Figure  2  for  a 
solid  aluminum  target.  Initially  most  of  the  energy  in  the  target  is  thermal.  As  the  laser  energy  is  deposited, 
the  thermal  and  ionization  energies  rise.  Kinetic  energy  increases  r£q)idly  as  the  ablation  process  becomes 
important,  and  is  comparable  with  the  thermal  energy  at  the  peak  of  the  laser  pulse.  Because  the  timescale 
is  so  short,  the  radiation  energy  is  about  three  orders  of  magnitude  less  than  the  thermal  energy  at  0.40 
picoseconds.  At  late  times,  the  plasma  expands,  cools  and  recombines.  Eventually,  much  of  the  ionization 
energy  will  be  radiated  away  (radiative  recombination).  Thus,  the  x-ray  radiation  pulse  from  the  target  will 
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be  much  longer  than  the  incident  laser  pulse.  The  energy  histories  for  the  layered  target  configurations 
closely  approximate  Figure  2. 


Figure  1.  IncidentlaaerpoweriKnmaUzedtothepeak 
of  the  laser  pulse  (lO^^IVycm^),  absorption  fraction, 
and  temperature  at  the  critical  surface  in  tens  of  keV 
are  shown  as  functions  of  time  for  P-wave  dqrasition. 


Figure  2.  Energy  history  Iot  a  10^^  W/cm’  pulse 
incident  on  a  solid  aluminum  target  The  total  energy 
in  the  plasrru  is  shown,  as  well  as  the  thermal,  kinetic, 
ionization  (and  excitation)  and  radiated  energies. 


Figure  3  shows  the  electric  and  magnetic  fields  and  the  local  energy  deposition  rate  in  the  target  near  tl^ 
peak  of  the  laser  pulse  (0.30  picoseconds)  for  an  aluminum  target.  The  laser  pulse  is  assumed  incident  from 
right  to  left  The  fields  are  averaged  over  an  optical  t^cle,  and  normalized  with  respect  to  the  amplitudes 
of  the  incident  electric  and  magnetic  fields.  The  magnetic  field  is  plotted  on  a  linear  scale.  However,  the 
electric  field  becomes  so  large  near  the  critical  surface,  that  the  logarithm  (logio)  of  the  field  is  plotted. 
The  location  of  the  critical  surface,  i.e.  where  the  incident  laser  frequency  is  equal  to  the  plasma  frequency 
(u  =  wp)  is  represented  as  Rent-  The  density  is  plotted  in  a  stepwise  manner  to  show  the  spatial  resolution 
of  the  hydro  mesh. 

In  the  absence  of  an  energetic  electron  beam,  the  temperature  profile  would  be  mostly  dominated  by 
thermal  conduction  from  the  laser  deposition  region  and  would  exhibit  a  rapid  decay  into  the  cold  aluminum. 
The  density  is  nearly  constant  behind  the  blowoff,  although  a  small  disturbance  (a  rarefaction  and  a  weak 
shock)  can  be  seen  near  x  =  1.0  in  Figure  4  due  to  the  fast  electron  heating.  Radiative  losses  are  the 
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greatest  near  x~2.0n  where  the  plasma  is  hot.  A  secondary  cooling  peak  occurs  at  about  0.90  m.  but  is  is 
not  in  the  range  of  the  plot. 


Figure  S.  Enagy  DqxMitioii  into  an  Aluminum 
Tnget  The  rale  of  eneigy  dqwdtion  into  the  taiget 
dSitldt  is  plotted,  and  the  locatimi  of  the  critical  sur- 
fve  Rerit  is  shown.  The  local  electric  and  magnetic 
helds,  averaged  over  an  (qrtical  cycle,  are  shown  nor¬ 
malized  to  the  anq>litude  of  the  incident  fidds,  E„  and 
Bo.  Mass  density  p  is  plotted  in  a  stq>wise  manner. 


Figure  4.  Hydrod]mamic  profiles  in  the  target  at  0  JO 
picoseconds.  Mass  density  p  [p/cm’],  temperature  T 
[eV\,  pressure  P  [erp«/cm’zlO~^’],  and  radiative 
cooling  dEjt/dt  [ergt/paee  —  cm’z  are  plotted 
as  functions  of  position. 


The  laser  (dEt./dt)md  fast  electron  (dEr/dt)  energy  deposition  profiles  are  represented  in  Figure  5. 
The  local  Ka  powers  by  ionization  state  produced  by  the  deposition  of  fast  electrons  are  also  ^own, 
including  neutral  (I)  through  boron-like  aluminum  (DQ.  Most  of  the  neutral  and  magnesium-like  aluminum 
Ka  radiation  originates  in  the  cold  target  region  lying  between  about  1.0  p  and  1.7  p.  Intermediate-ion  Ka 
radiation  (firom  fluorine-like  through  sodium-like  aluminum)  originates  in  the  hot  laser-deposition  region 
(x  fa  1.8  p)  and  in  the  warm  fast-electron  heated  region  (x  »  0.9  p).  The  Ka  radiation  from  highly  ionized 
aluminum  (oxygen-like  aluminum  and  above)  comes  mostly  from  the  hot  laser  deposition  region. 

The  instantaneous  distribution  of  Ka  radiation  is  correlated  with  the  local  atomic  populations.  At  early 
times  Giefore  significant  ionization  can  take  place)  only  Ka-I  and  Ka-II  is  seen.  Most  of  the  Ka  firom  highly 
ionized  aluminum  is  created  at  intermediate  times,  near  and  shortly  after  the  peak  of  the  laser  pulse,  when 
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tlw  plasma  reaches  maximum  ionization.  At  later  times,  the  plasma  cools  and  recombines.  Even  spatially 
integrated  Ka  radiation  can  give  a  wealth  of  information  about  the  instantaneous  state  of  the  plasma. 

The  instantaneous  spatially  integrated  x-ray  emission  spectrum  for  the  K-shell  and  K<,  lines  is  plotted 
as  a  function  of  photon  energy  at  0.30  picoseconds  (near  the  peak  of  the  laser  pulse)  in  Figure  6.  The 
spectral  lines  are  identified  on  the  figure.  The  most  distinguishing  feature  of  the  spectrum  is  the  series  of 
lines  generated  by  the  energetic  electrons.  The  lines  from  the  four  lowest  aluminum  ionization  stages 
(Ka-I  through  Ka-IV)  lie  so  close  together  in  energy  that  they  appear  in  the  spectrum  as  a  single  line;  the 
remaining  Ka  lines  ate  well  resolved.  The  Ka  I-TV  line  is  prominent  in  the  spectnun,  and  is  comparable  in 
intensity  with  the  He-a  line. 


Dlipltconeat  (p) 


Figure  5.  Enogy  dq)osition  at  0.30  picoseconds 
for  P-polaiization.  The  logarithm  of  the  rate  of  laso' 
eneigy  dqwsition  dEL/dt  and  fast  electron  eneigy  de¬ 
position  dEr/dt  [xlO^^’istF/cm^]  are  plotted  with  the 
Ka  powers  for  various  ionized  states  of  aluminum  pro¬ 
duced  by  the  dqwsition  of  fast  electrons.  Th^  include 
neutral  (I)  through  boron-like  aluminum  (DQ. 


lBcr|7  (ktf) 


Figure  0.  Emission  spectrum  at  0.30  picoseconds 
for  an  aluminum  target,  showing  magnitudes  of  the  Ka 
and  K-sbell  lines.  The  He-o  line  is  comparable  to  the 
neutral  Ka  line  of  aluminum. 


The  relative  intensities  of  the  lines  reflect  the  local  atomic  populations  which  the  energetic  electrons 
encounter.  As  Figure  S  illustrates,  high  states  of  ionization  are  encountered  in  the  blowoff  plasma,  and 
low  ionization  states  ate  encountered  elsewhere.  A1  n  and  A1  HI  are  the  most  abundant  ion  species  in  the 
cold  target,  and  neutral  aluminum  is  down  in  abundance  by  about  an  order  of  magnitude.  Most  of  the  fast 
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electrons  are  generated  near  the  critical  surface;  since  the  stopping  power  varies  with  electron  energy  as 
ln(E)/E,  K-shell  vacancies  are  created  more  efficiently  about  a  micron  into  the  target,  when  the  electrons 
have  slowed  down  (A1  n  through  A1 V  dominates  in  this  region).  Finally,  A1  DC  dominates  the  hot  blowoff 
region;  the  lower  ionization  states  are  substantially  depleted.  In  the  simulations,  Ko  transitions  for  A1  X 
and  higher  are  not  considered.  The  result  of  these  considerations  is  seen  in  the  relative  intensities  of  the  Ko 
lines  in  Figure  6. 

3^  Layered  Silicon/Aluminum  Target 

While  the  relative  intensities  of  the  lines  can  provide  detailed  information  about  atomic  populations, 
it  is  difficult  to  use  this  knowledge  in  the  differentiation  of  fast  electron  phenomena  from  laser  absorption 
effects.  However,  the  use  of  multiple  materials  in  the  target  with  layers  of  varying  thicknesses  provides 
additional  spectral  information.  The  range  of  energetic  electrons  in  the  target  can  be  determined  directly 
from  such  experiments,  and  the  electron  energy  distribution  function  can  thus  be  inferred.  In  addition,  it 
may  be  possible  to  make  quantitative  measurements  of  laser  absorption  as  a  function  cf  position  in  the 
target,  and  to  verify  assumptions  regarding  multiphoton  processes,  etc. 

A  series  of  layered  target  simulations  was  carried  out,  where  the  initial  conditions  were  essentially 
identical  with  the  aluminum  slab  target  described  above  (0.22  psec  FWHM  duration  Gaussian  pulse  incident 
at  30  degrees  from  the  normal  with  peak  power  1.0x10”  W/cm*  and  a  wavelength  of  0.26  /*)•  In  the  first 
case,  the  laser  was  incident  on  a  0.10  silicon  layer  coating  a  solid  density  aluminum  substrate.  The 
energetics  of  the  deposition,  the  energy  history  and  the  hydrodynamic  response  of  the  target,  and  the 
production  and  deposition  of  energetic  electrons  in  the  target  are  approximately  the  same  as  for  the  case  of 
the  plain  aluminum  target,  as  shown  in  Figures  1  through  5.  Naturally,  there  are  differences  in  the  details, 
because  silicon  has  a  different  solid  density  and  different  ionization  potentials,  etc.  However,  because 
aluminum  and  silicon  are  adjacent  elements  on  the  periodic  table,  these  differences  are  not  large.  The  most 
significant  difference  between  the  simulations  is  uie  resulting  radiation,  particularly  the  K-shell  emission. 

Figure  7  shows  the  K-shell  emission  spectrum  at  the  peak  of  the  laser  pulse  for  the  case  of  a  0.10 
silicon  coating  on  an  aluminum  target.  As  in  the  plain  aluminum  case  discussed  above,  the  fast  electrons 
are  assumed  to  be  monoenergetic  at  10  keV.  Now  there  are  two  sets  of  Ka  radiation  in  evidence.  Clearly, 
the  aluminum  plasma  is  not  as  hot  as  in  Figure  6;  aluminum  K-shell  radiation  is  completely  absent  except 
for  the  Ka  lines,  and  the  K*  V  line  is  almost  two  orders  of  magnitude  down  from  the  I-IV  line.  However, 
the  silicon  radiation  suggests  that  the  outermost  0.10  micron  of  plasma  is  hot.  Some  care  must  be  taken 
in  the  identification  of  the  lines,  because  certain  aluminum  lines  are  almost  coincident  with  silicon  lines. 
The  Rydberg  series  lines  from  H-  and  He-like  silicon  are  observed;  these  result  from  the  decay  of  np  states 
(n  <  5)  to  the  Is  orbitals.  The  silicon  Ka  radiation  q)parently  comes  from  two  distinct  "ources.  The  Ka  I-V 
line  originates  in  cooler  plasma  (behind  the  critical  surface),  and  the  Vn  through  X  lines  come  from  the 
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hot  region  just  ahead  of  the  critical  surface.  These  regions  correspond  to  those  described  in  Figure  5.  At 
this  time,  the  back  edge  of  the  silicon  layer  just  barely  extends  into  the  cool  region  extending  from  about 
1.2  /i  to  1.8 

If  the  fast  electrons  have  less  energy,  their  range  is  correspondingly  decreased,  and  the  plasma  just 
behind  the  critical  surface  becomes  hotter  (the  electron  energy  is  deposited  in  a  smaller  volume).  Figure 
8  shows  the  K-shell  emission  spectrum  for  an  identical  simulation,  except  that  the  fast  electron  energy  is 
assumed  to  be  4.0  keV.  The  aluminum  still  does  not  get  hot  enough  to  produce  thermal  radiation  in  the 
K-shell;  however,  the  signature  is  considerably  different.  The  Ka  V  and  VI  lines  are  nearly  as  strong  as 
the  I-IV  line,  and  the  Vn  line  is  almost  an  order  of  magnitude  less  intense.  The  I-IV  line  is  substantially 
less  intense  than  in  Figure  7.  The  silicon  emission  is  also  altered,  but  in  a  more  subtle  way.  All  of  the 
silicon  Ka  lines  are  more  intense,  because  more  fast  electrons  are  being  stopped  in  the  silicon  layer.  In 
addition,  the  reiat^ve  brightness  of  the  Ka  VI  line  is  increased,  since  the  plasma  immediately  behind  the 
critical  surface  is  wanner. 

The  local  Ka  powers  by  ionization  state  produced  by  the  deposition  of  fast  electrons  for  the  present 
case  is  plotted  as  a  function  of  position  in  the  target  in  Figure  9.  The  portion  of  the  target  from  1.2  to  2.2 
microns  is  shown,  since  the  fast  electrons  have  substantially  reduced  range  at  4.0  keV.  The  plasma  at  and 
ahead  of  the  critical  surface  looks  qualitatively  similar  to  the  corresponding  region  in  Figure  5.  The  details 
are  different  because  the  fast  electron  deposition  is  altered  and  because  the  material  is  silicon  rather  than 
aluminum.  There  is  substantial  production  of  silicon  Ka  DC  and  X  in  the  hot  plasma.  However,  deeper  into 
the  target  (between  about  1.6  and  1.8  microns),  the  plasma  is  qualitatively  different.  It  is  hotter  because 
of  the  reduced  energetic  electron  range,  and  the  ionization  state  is  higher.  The  fast  electrons  are  creating 
K-shell  vacancies  in  more  highly  ionized  material  (predominately  aluminum).  There  are  strong  gradients 
in  the  Ka  powers  near  1.78  microns,  which  is  at  the  aluminum/silicon  interface.  The  plot  illustrates  that 
the  actual  Ka  production  is  more  complicated  than  the  spectra  suggest.  There  are  at  least  six  regions  in  the 
target  plasma  where  the  ionization  physics  is  distinct.  Deep  in  the  target,  the  (aluminum)  plasma  is  cold 
and  essentially  unionized.  Between  about  1.6  and  1.78  microns  the  (aluminum)  plasma  is  warm  (of  the 
order  100  eV),  and  aluminum  II-VI  predominates.  Sandwiched  between  this  and  the  critical  surface  region, 
the  plasma  is  also  of  the  order  100  eV,  and  silicon  IV- VI  is  most  prevalent.  Near  the  critical  surface,  the 
silicon  is  hot  and  dense;  silicon  Vm-X  predominates.  Ahead  of  this  region  (near  1.9  microns)  the  density 
is  somewhat  lower  and  silicon  IX-X  is  most  abundant.  Since  we  only  consider  Ka  lines  of  silicon  through 
boron-like,  the  plot  contains  no  information  about  silicon  XI-XV.  Finally,  near  2.2  microns,  intermediate 
ionization  states  of  silicon  dominate,  since  the  density  is  low  and  *he  plasma  is  far  from  equilibrium.  In 
this  way,  the  fast  electron  energy  spectrum  can  be  estimated.  The  layer  thicknesses  can  be  varied  to  give 
additional  information  about  the  range  of  the  electrons.  The  results  of  such  simulations  will  be  incorporated 
into  a  subsequent  paper. 
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Figure  7.  Emission  spectnim  at  0.30  picoseconds  Figure  8.  Emission  spectrum  at  030  incoseconds 

for  an  aluminum  target  with  an  0.10  p  alicon  coating.  for  an  aluminum  target  with  an  0.10  p  silicon  coating. 

The  K-sbell  (including  the  Ka  lines)  of  aluminum  and  The  fast  electrons  are  assumed  to  be  monoenergedc  at 

silicon  are  shown.  The  fast  electrons  are  assumed  to  4  keV. 

be  monoeneigetic  at  10  keV. 

The  use  of  this  diagnostic  is  further  illustrated  by  conqiaring  the  spectra  from  aluminum/silicon  and 
silicon/aluminum  taigets.  Figure  10  shows  the  K-shell  emission  for  the  case  of  a  0.10  p  aluminum  coating 
on  a  silicon  substrate.  The  laser  characteristics  are  the  same  as  for  the  previous  cases,  and  the  energetic 
electrons  are  taken  to  be  at  4.0  keV.  Now  aluminum  K-shell  radiation  (the  Rydberg  series  lines  from  H- 
and  He-like  aluminum)  is  in  evidence,  and  silicon  K-shell  radiation  is  absent  except  for  the  Ka  lines.  The 
aluminum  H-a  line  is  situated  just  below  silicon  Ka  sequence  in  energy  (and  can  easily  be  mistaken  for  the 
first  line  of  the  sequence);  it  is  at  least  an  order  of  magnitude  more  intense  than  the  silicon  Ka  I-V  line.  The 
Ka  lines  from  higher  stages  of  ionization  in  silicon  (VI  and  VII)  fall  off  rq)idly.  This  result  is  in  contrast 
with  the  silicon-on-aluminum  spectrum,  where  the  first  three  Ka  lines  were  comparable  in  intensity.  This 
implies  that  the  silicon-on-aluminum  plasma  is  slightly  warmer  behind  the  aluminum/silicon  interface. 
This  can  be  explained  by  the  different  energetic  electron  stopping  powers  of  aluminum  and  silicon  and  by 
the  different  electrical  conductivities  of  the  species  which  can  produce  variations  in  the  laser  deposition. 

The  aluminum  Ka  lines  are  all  comparable  in  intensity;  however,  the  Ka  Vm  line  is  weak.  As  before, 
this  implies  that  the  radiation  is  coming  from  two  distinct  regions.  The  Ka  I-Vn  is  coming  from  a  warm 
region  just  behind  the  critical  surface.  The  Ka  DC  comes  mainly  from  the  hot  plasma  at  and  ahead  of  the 
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critical  surface.  The  spectra  in  Figures  8  and  10  thus  convey  similar  information.  The  differences  between 
them  illustrate  the  complications  that  arise  in  interpreting  spectral  data. 


PitpUcMtat  (|i) 


btrfj  (ktV) 


Figure  9.  Ka  powers  for  various  ionized  states  of 
aluminum  and  silicon  produced  by  the  dqposition  of 
fast  dections.  They  include  neutral  (I)  through  boron¬ 
like  aluminum  (DQ,  and  neutral  through  boron-like 
silicon  (X). 


Figure  10.  Emission  spectnun  at  0.30  picoseconds 
for  a  silicon  target  with  an  0.10  n  aluminum  coating. 
The  fast  elections  are  assumed  to  be  monoeneigetic  at 
4keV. 


4.  CONCLUSIONS 

The  purpose  of  this  investigation  was  to  simulate  the  radiation  hydrodynamics  of  an  ultrashort  pulse 
laser  produced  plasma  and  to  study  how  the  emitted  x-rays  could  provide  information  about  the  energy 
(^position  in  the  target,  the  production  and  transport  of  energetic  electrons,  and  the  resulting  ionization  of 
the  plasma.  The  model  which  was  used  included  a  relatively  detailed  treatment  of  the  atomic  populations 
and  radiation.  The  laser  deposition  was  performed  in  a  self-consistent  maimer  by  solving  the  appropriate 
Helmholtz  equation.  The  treatment  of  the  energetic  electrons  incorporated  a  number  of  prescriptions 
(production  rate  and  energy  distribution),  but  it  provided  for  the  self  consistent  transport  of  the  electrons 
and  the  production  of  K-shell  vacancies.  One  aspect  of  the  calculation  that  needs  improvement  is  the 
spatial  resolution  in  the  vicinity  of  the  critical  surftK».  Since  Ka  radiation  is  produced  almost  entirely  by 
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the  energetic  electrons,  these  inner  shell  K-lines  are  particularly  useful  in  providing  information  on  energy 
transport  due  to  fast  electrons.  However,  the  ratios  of  the  Ka  lines  reflect  the  local  atomic  populations 
where  the  electrons  produce  K-shell  vacancies.  It  is  challenging  to  accurately  determine  the  electron  energy 
spectrum  from  this  data  alone.  The  details  of  the  laser  deposition  and  other  factors  which  can  alter  the 
density  and  tempoature  profiles  (and  the  local  ionization  of  the  plasma)  can  also  affect  the  line  ratios. 
The  deposition  is  controlled  by  the  local  electrical  conductivity,  and  progress  continues  to  be  made  by  a 
number  of  researchers  in  developing  a  satisfactory  conductivity  naodel  for  the  conditions  encountered  in 
these  plasmas.  The  use  of  very  thin  multiple  material  layers  (as  has  been  explored  by  Chm  et  al.  1993) 
is  one  way  to  provide  additional  information  about  the  ionization  of  the  plasma  as  a  function  of  position 
in  the  target,  and  thus  distinguish  between  these  factors.  It  is  hoped  that  the  analysis  of  the  dynamics  of 
layered  targets  contained  in  the  last  section,  although  far  from  exhaustive,  proves  useful  to  experimentalists 
in  interpreting  their  data  We  also  investigated  tl^  possibility  of  plasma  conditions  being  generated  in 
the  target  that  are  conducive  to  the  creation  of  population  inversions.  For  relatively  high  laser  intensities 
(S-wave  deposition  at  3xlO^'W/em’,  for  exaiiq)le),  ten^perature  and  density  profiles  in  the  blowoff  plasma 
have  tentatively  been  found  to  support  inversions  in  a  number  of  upper  states  of  lithium-like  aluminum. 
The  most  favorable  conditions  seem  to  occur  for  electron  densities  of  the  order  10^em~’.  These  results 
will  be  presented  and  discussed  elsewhere. 
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Short-pulse  Driven  Selenium 
X-ray  Lasers 


K.  G.  Whitney,  A.  Dasgupta,  and  P.  E.  Pulsifer 


Plasma  Physics  Division 


Abstract:  Ihe  of  intense  picosecond  heating  and  time  depeaSeas.  ionization  on  the 

population  dtynamics  of  Ne4ikB  sdenhim  is  investigated.  We  find  that  highly  non-eqoilibdnm 
fdasma  conditions  and  a  laige  gain  ooefiKdent  >  600  cm~^  can  be  generated  for  the  JM)-1  line 
at  182  k  at  KiF  laser  intenaties  >  3  x  10**  W/cm*.  Thns,  die  observation  of  x-ray  lasing  at  or 
near  the  critical  sarfaoe  of  a  KrFdiort-polse,  laser  driven  [dasmagycars  feasible.  Its  observation 
would  provide  important  tests  of  dieotetical  modds  im  laser  absoiption,  transient  ionization,  and 
x-n^  laser  pnmping. 
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Bariy  in  tiwir  development,  it  was  conjectured  that  high  power,  picosecond  pulsed  lasers 
could  become  useful  table>tt^  drivers  for  x>ray  lasers.^  A  number  of  different  schemes  for 
generating  population  inversions  in  die  plasmas  heated  by  picosecond  or  subpicosecond  lasers 
were  then  prqposed.^*’  In  dus  letter,  we  present  the  results  of  calculations  that  demmistrate  die 
feasibility  of  generatiiig  large  transient  populadrm  invetsimis  widun  die  neon-like  imuzadon  stage 
of  selenium  using  diott-^iulse  picosecond  lasers.  These  inversimis  are  large  enough  to  produce 
a  sidnrated  amplified  qiontaiieous  emission  (ASE)  firom  a  plasma  <  1  mm  in  length.  These 
x-ray  laser  experiments  could  be  carried  out  with  1  to  10  picosecond  KrF  laser  pulses  heating 
a  pve-fotmed  selenium  plMma  at  or  near  the  ctidcal  surface.  The  use  of  selenium  in  these 
eaqieriments  has  two  advantages.  One,  neon-like  selenium  x-ray  lasers  have  already  been  operated 
in  long-imlse  (~  500  ps)  sr*'^****^  plasmas,^  and  a  number  of  dieir  properties  have  been  wdl 
studied  eaqienmentally.***  TWo,  fairly  oompteheiisive  models  of  their  population  tfynamics  have 
been  oonstracted  and  benchmarked  against  esqierimental  observations;^'*  one*  isusedinobtaitiing 
the  calcnlational  results  described  in  dus  letter. 

Shoct-pulse  generated  plasmas  ate  also  of  general  interest  as  picosecond-pulsed  sources  of 
kflovolt  z  rays.*  Bqrerimental  obsetvatimis,  in  dus  case,  dearly  indicate  die  need  totrapidly 
heat  a  pre-foimed  plasma  in  order  to  generate  die  maximum  rmmbef.ofznqrs.  The  calculations 
we  present  support  dus  observation.  The  faydrodtynamiceaqMuision  of  Uowoffplasinas  generally 
occurs  on  die  order  of  mmosecond  time  scales;  thus,  die  picosecond  heating  and  ionization  of 
a  pre-fpnned  itrignimn  plama  can  initially  be  studied  by  ignming  the  effects  of  hydrotlynamic 
egpanann-  For  dus  purpose,  we  the  hot-spot  modd  described  eadier^*  to  determine 

die  inwiyjifinn  tesponse  of  arfemiim  to  picosecond  heating.  The  sdenium  atomic  modd  in  this 
y!^icni5irio«i  consists  of  die  faydrogenic  lamped  state  modd  described  in  Ref.  (11),  to  vduch  the 
neon-like  and  fluorine-like  structure,  wfaidi  is  described  in  Ref.  (8),  was  added.  In  this  model, 
a  full  descn^m  is  given  of  the  transient  exdtadon  and  ionizatitm  of  all  26  r¥=3  states  in  die 
neon-like  imiizadon  stage  from  which  x-ny  lasing  has  been  observed.  Specifically,  lasing  was 
seen  in  two  J=0-1  and  two  J=2-l  3p  to  3s  transitions.  However,  at  the  plasma  density  of  our 
problem,  which  will  correspond  roughly  to  the  ion  density  of  a  KrF  laser  critical  surface,  the  J=2-l 
transitions  cannot  be  inverted  and  show  no  gain  even  in  transient  ionization.  The  J=0,  3p  states, 
on  the  other  hanH,  are  strongly  populated  from  both  the  ground  state  and  from  the  J=l,  3d  states 
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due  to  the  huge  monopole  and  dipole  exdtatitm  rates  from  the  ground  state  to  both  the  J=0, 3p  and 
Jal,  3d  states,  respectively.  Hence,  during  transient  excitation,  a  large  populatitm  invertion  can  be 
estaWudted  between  the  JsO,  3p  and  Jsl,  3s  states  at  plasma  (tensities  correspcMiding  to  the  KrF 
critical  surface. 

In  die  hot-qxM  mcxtel,  the  following  ntmlineariy  coupled  equations  are  solved: 

^  =  (* = 1 . ",  (1) 


M  M 

+  Ei^=S,-Qk-Ri- -  ^0’ 

idH-SvC'-'). 


(2) 

(3) 

(4) 


vdiere5c  rqnesents  die  laser  energy  absoiptionratepernnit  volume,  and  isajdienomenological 

iqptesentation  of  die  heat  conductitm  losses  dial  mighi  occur  from  die  d^  focus  of  a  laser-heated 
plasma.  The  seleniiim  model  contams  iisl89  states  and  popolaricm  denstdes,  {iV)*}.  {Z^}  and 
{E^}  are  die  icmic  (diarges  and  energies  of  these  states.  Thejr  are  coupled  by  the  cmnplefe  set 
(tf  collisioiial  ionizadtm  and  exdtadon  rates  and  ladiadve  and  coUisicmal  decqr  rates, 
described  in  Ref.  (8).  The  equadon  of  modon,  dEifdt  =  Qi  —  iZ,,  is  derived  from  Eqs.  (1).  Ri 
is  die  sum  of  opdcally  diin  line  and  free-bound  radiadon  loss  rates,  which  are  obtaiiied  from  the 
radiation  terms  in  die  rate  eqoad(ms(E(]s.(l)).  Q{  is  die  rate  ofinelasdccollisi(»al  energy  transfer 
from  the  electrons  to  the  ions,  vduch  is  balanced  in  collisitmal-radiative  ecjuilibrium  (CRE)  byRi. 
Rhrem»  ts  die  rate  of  bremsstrahlung  losses  from  the  plasma.  In  picosecond-heated  plasmas,  both 
Ri  and  R^remt  are  negligble  relative  to  5*,  Qh,  and  the  rate  of  thermal  oiergy  transfer  from  die 
dectrons  to  die  ions,  given  by  die  last  term  cm  the  right  hand  side  of  Eq.  (4),  which  is  proportional 
to  Te  -  Ti.  The  expressicm  for  diis  teim  is  taken  frmn  Braginskii.^^  It  is  inversely  proportional  to 


the  dectnm-ion  collision  time  . 

The  expressitm  for  taken  from  Ref.  (10),  like  5e,  is  proportional  to  Qh  is  also 
determined  in  terms  of  a  free  parameter,  the  dectron  temperature  gradient  scale  length,  ig.  This 
scde  length  allows  one  to  phenomenologically  estimate  the  maximum  effect  of  heat  flow  energy 
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losses  <xi  tile  time  Ustoiy  of  the  election  temperature,  Te.  These  losses  are  averaged  over  the 
volnme  of  an  absorption  qiheie,  i.e.,  a  ^tiiere  ^(diose  radius,  r,  is  rou^y  of  tiie  order  of  the 
inverse  biremsstrahlnng  absmption  length  of  the  laser  radiation  near  tiie  critical  surface.  Thus,  tiie 
expression  ftn  used  in  tiiese  calculations  is: 


Qk  = 


3K«r« 

rf. 


m, 


STe 

Tie’ 


(5) 


udieneverQ^  is  less  tiumafiee  Streaming  limit,  =  O.ZNekTeVtkm*y9beKvthm  =  y/2kTe/mt. 

Otiieratise,  Qa  is  set  equal  to  Q{“. 

We  consider  tiie  {di^cal  situation  in  vtiudi  a  premise  is  used  to  generate  a  blowoff 
jtiasma  containing  i<m  densities  tiut  will  produce  a  critical  surface  coneqionding  to  tiie  0J2S  /tm 
wavelengtii  of  a  Kri’  laser  when  fhDy  ionized.  Hie  electron  density  of  a  0.2S  foa  wavelengtii 
laser  is  =  1.8  x  10”.  The  main  KiP  pulse  can  then  be  used  to  resonantly  heat  and  ionize 
sdeninm  near  this  surface  where  tiie  ion  densi^  is  <  7  x  10^".  At  tins  density,  tiie  gains  of 
tiie  J=0-1  and  Ja2-1  lines  at  1  keV  would  be  eitiier  zero  or  near  zero,  if  the  {tiasma  were  in 
(CRE).*  We  assume  tiiat  tile  blowoffplasma  has  a  suffidentiy  small  scale  lengtiitiiat  a  substantial 
portion  of  tiie  main  laser  pulse  will  tiien  be  absorbed  at  or  near  tiie  cntical  surface  Ity  resonant 
inverse  bremsstralilniig.^’  Near  tiie  critical  surface,  laser  absorption  occurs  at  a  rate  given  by 
5«(i)  rs  ic£(f)/i;(i),  ediere  /l  is  tiie  inddent  laser  intensity  and  ki;  is  the  absorption  coeffident 
given,  fitom  Refs.  (10)  and  (13),  by 


1  8.3  X  10»i;f^*(l  - 

KL  ~  (iV2/N,)lnA(JVeW) 


(6) 


In  A  is  the  Coulomb  logaritm  and  AT*  is  tiie  critical  electron  density  corre^nding  to  0.25  paa 
radiation. 

(l)-(6)  describe  tiie  resonant  inverse  Bremsstrahlung  heating  of  a  small  volume  dement 
of  plasma  whose  dectrtms  ate  cooled  by  botii  dectnm-ion  dastic  and  indastic  collisitxis  and, 
itiienomenologically,  by  (potentially  large)  flux  limited,  heat  conduction  energy  losses.  The  volume 
dement  used  in  our  calculations  has  a  fixed  itm  density  of  6  x  10^**  and  a  fixed  radius  of  10  /xm, 
^dnch  is  of  the  order  of  the  laser  absorption  lengtii.  This  volume  dement  is  heated  with  a  KrF 
laser  pulse  of  2  picoseconds  FWHM.  Figs.  (1)  and  (2)  show  a  typical  response  of  the  temperature 
and  ionization,  reflectively,  to  a  Gaussian  laser  pulse  of  peak  intensity,  =2x10^®  W/cm^ . 
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In  (»der  to  illustrate  tlw  nq}idity  of  the  ionizadoo,  we  began  these  calculations  with  a  cold  (~  2 
eV)  weakly  itmized  plasma.  The  laser  pulse  peaks  at  3  pcoseconds,  and  the  scale  loigdi  for  heat 
conduction  losses  was  initially  taken  to  be  6  k  so  that  heat  losses  were  flux  limited  for  roughly  the 
first  two  picosectmds  of  the  calculaticm. 

Because  large  heat  losses  are  simulated  in  this  calculatimi,  the  electron  temperature,  which 
had  apeak  above  1400  eV  at  4  ps,  fiell  qniddiy  to  700  eV  at  lOps.  The  rise  in  die  ion  temperature 
as  wdl  as  die  iqiid  ionization,  isduch  reached  die  neon-like  ionization  stage  in  6  ps  (Hg.  (2)), 
also  contributed  to  die  cooling  of  the  electrons,  but  to  a  lessor  degree.  However,  in  spite  of  diis 
oodUng,  die  idasma  continued  to  ionize  over  die  entire  10  ps  interval  because  the  rise  in  ionization 
lagged  die  rise  in  T*  and  because  r«  readied  such  a  large  peak.  Note  diat,  at  6  x  10^°  ions/cm^, 
adenhun  Ionizes  nqddly  to  die  neon-like  ionization  stage  before  the  ionization  is  slowed  by  die 
much  larger  exdtation  energies  of  die  L-dieU  compared  to  die  M-shelL  At  6  x  10^°  ions/cm’, 
dieiefote,  selemum  is  an  ideal  candidate  material  for  generating  population  inversioos  widiin  die 
NoJike  ionization  stage  on  the  10  iMcoseoond  time  scale  of  this  calculation. 

Because  the  ion  temperature  is  low,Dop|derbn>adeningis  only  asmallconectiwi  to  ctdlisional 
broadening  of  die  laang  lines  in  diese  calculations.  Therefoite,  ^^aigt  profiles  were  needed  to 
compute  the  gains  of  die  JssO-l  lines  at  line  center^*  s  (ilA^/8ar)(iV«  —  {gM/gi)Ni)IAt/Vt 
where  A  is  the  radiative  decqr  rate  for  die  transition,  A  is  die  radiated  wavelengdi,  and 
{Ni^gt}  are  the  populations  and  d^enera^  factors  of  die  t^per  and  lower  states  of 
die  transition  reflectively.  Frv  die  value  of  a  Yoigt  profifo  at  line  center,  l/Auv^  we  used  an 
expression  derived  by  J.  .^iruzese:^* 


1  ^  l+g _ }_ 

Ai/v  1  +  1.66o  -I-  1.772o*  ^/5rAi/u  ’ 


(7) 


where  ud  =  V2kTi/{mic*){c/X)  is  the  Doppler  linewiddi  and  a  =  Tt/{^^Aud)  is  the  Voigt 
parameter.  Fr  is  die  sum  of  die  inverse  lifetimes  (radiative  plus  colliaonal)  of  the  upper  and 
lower  states.  In  our  atomic  model,  these  linewiddis  are  fiproximately  100  mA  for  both  J=0-1 
lines.  Moreov^,  die  calculated  J=2-l  linewiddis  were  roughly  lA  times  larger  than  diose  of  the 
J=0-1  lines.  These  results  benchmark  reasonably  well  against  recent  experiments,'  in  which  a 
SO  mil  linewidth  was  measured  for  the  J=2-l  line  at  206  'A.  Furthermore,  it  was  speculated  that 
collisitmal  effects  were  playing  a  significant  role  in  determining  the  line  profile. 
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In  the  calculation  under  dicussion  (see  Figs.  (1)  and  (2)),  the  gains  of  the  J=0-1  lines  had 
rou^y  the  time  behavior  of  T.;  namely,  diey  acquired  nonzero  values  after  4  ps,  aiul  they  peaked 
at  S  J  ps.  Then,  they  steadily  fell  to  roughly  1/3-rd  of  their  peak  values  at  10  ps.  As  die  laser 
intensity,  /£**”,  of  the  2  ps  pulses  was  varied,  the  size  of  the  J=0-1  population  inversions  varied, 
depending  diiecdy  <m  die  amount  of  laser  energy  deposited  in  the  plasma.  This  behavior  is  shown 
mHg.(3).  Fori^^”  =2  x  10^‘,  corresponding  to  Figs.  (1)  and  (2),  the  maximum  gain,  of 
die  182  it  J=0-1  line  was  224  ,  sdnle  die  JsO-1  line  at  169  >1  adueved  a  maximum  gain  of  87 

cm~^.  However,  jrma«  continued  to  rise  as  was  increased.  Beyond  an  intensity  of  3  x  10^* 
W/cm^,  the  182  k  line  is  {Hedicted  to  teach  gains  in  excess  of 600  cm~^.  The  FWHM  of  the  ASE 
pulses  that  would  be  generated  in  diese  plasmas  is  predicted  to  be  2  ps.  The  propagation  distance 
of  a  2  ps  x-ray  pulse  widun  an  anqditying  medium  is  0.6  mm. 

One  proUem  associated  with  die  gain  modds  diat  have  been  constructed  to  date  to  describe 
the  ptqiolatitm  dynamics  of  die  sdeninm  neon-like  itmization  stage  is  that  they  predict  the  largest 
gain  for  die  JsO-1  line  at  182  k.  In  long  pulse  x-rity  laser  eiqietiments,  however,  larger  gains  were 
measured  for  die  Ja^-l  lines.  This  problem  was  discussed  in  Ref.  (8),  udiereit  was  noted  that  a 
factor  of  two  decrease  in  die  magnitnde  of  the  calculated  monopde  excitation  cross  sectioi^  which 
coiqiles  die  ne<m-like  ground  state  to  die  (upper)  3p  state  of  this  182  k  transitimi,  was  sufElcient 
to  reverse  the  relative  magnitudes  of  the  JsO-1  and  Js2-1  gains  in  CRE  at  an  ion  density  of  10^’. 
The  effect  that  such  a  cross  section  reduction  has  in  diese  short-pulse  calculations  is  riiown  by 
the  dashed  curve  in  Hg.  (3).  TUs  reduction  in  the  monopole  rate  produced  less  dian  a  factor 
of  2  reduction  in  the  peak  gain  of  the  182  k  line  because  this  line  is  strongly  pumped  from  die 
uppermost  oftheJsl,  3d  states.  Moreover,  dus  rate  rednctitm  did  not  significandy  change  die  gain 
of  die  169  k  line. 

Three  factors  combine  to  produce  die  large  gains  seen  in  Hg.  (3).  One,  die  calculations  are 
carried  out  at  an  im  density  that  is  larger  than,  Ity  over  an  order  of  magnitude,  the  densities  of  the 
gain  medium  inferred  in  the  Img-pulse  x-ray  laser  experiments.  Thus,  collisitmal  exdtatitm  rates 
are  also  an  or^  of  magnitude  larger.  T\vo,  the  plasma  ionizes  through  a  successitm  of  highly 
non-equilibrium  states.  Hence,  Ity  the  time  neon-like  ground  states  are  plentiful,  the  dectrons 
are  already  hot  enough  to  ionize  the  plasma  into  die  higher  ionization  stages  of  the  L-shell,  well 
beyond  the  neon-like  stage.  Under  these  conditions,  unlike  those  foimd  in  CRE  plasmas,  there  is  no 


45 


exponential  temperature  dependence  to  the  ground  state  excitation  rates,  and,  consequently,  they 
are  larger  than  they  would  be  in  CRE.  Finally,  the  pofmlations  of  the  3d  states  that  feed  the  J=0, 3p 
states  are  siqtpressed  in  CRE  relative  the  odier  metastable  3d  states  due  to  strong  radiative  decays 
to  tbffi  ground  state.  The  ojqxmte  is  true  during  transient  excitation,  i.e.,  strong  radiative  decays 
imply  that  the  collisional  excitation  rates  of  these  3d  states  are  also  ^rtmg.  (In  our  calculations, 
tb^  ate  roughly  tme  to  two  orders  of  magnitude  laiger  than  die  excitadcm  rates  to  the  odier  ten, 
metastaUe  3d  states.)  Thus,dietwo3dstatesdiatstrongly  cotqilecoUisimiaUy  totbeSpJsOstates 
are  relatively  oveipt^Milated,  rather  than  undeipc^iulated,  and  dieir  role  in  populating  the  3p  states 
is  significandy  enhanced  compared  to  CRE  calculatitxis. 

Hgs.  (1H3)  define  the  intenaty  and  polsewiddi  requirements  for  a  short-pulse  laser  diat 
would  be  needed  to  cany  out  die  above  x-n^  laser  eiqienments.  Fig.  (4)  roughly  defines  the 
energy  requirements  of  the  laser.  It  shows  the  ammmt  of  energy  per  unit  mass  that  was  absorbed 
in  each  of  the  above  calculaticms.  For  examine,  a  cylinder  10  Atm  in  radius  and  05  mm  in  length 
containing  6  x  10^^  ions/cm^  contains  1.25  x  10~^/tg  of  selenium.  According  to  Hg.  (4),  heating 
dus  cylinder  widi  a  2  ps  laser  pulse  at  3  x  10^*  W/cm^  wouldreqniteS.4J  of  energy.  In  addition 
to  dus  energy,  one  must  factor  in  die  laser  energy  absorbed  in  odier  portions  of  the  blowoff^lasma 
as  well  as  reflected  energy. 

At  the  laser  intensities  used  in  these  calculatitms,  highly  non-Maxwellian  elecmm  distributions 
will  be  ynerated  by  intensi^  of  die  inverse  Bremsstrahlung  heating.^*'^^*^*  At  ai^  given  level  of 
ionization,  diese  distribations  will  enhance  the  rates  for  exciting  the  plasma,  and,  thus,  th^  will 
lead  to  larger  pqpulaticm  inversioiis  than  were  attained  for  die  calculated  gains  shown  in  Fig.  (3). 
Hence,  diere  is  reastm  to  believe  that  the  calculated  gains  shown  in  Fig.  (3)  are  smaller  than  may 
be  achieved  in  an  actual  short-pulse  KrF  experiment  Anodier  effect  of  die  non-Maxwellians  will 
be  to  alter  the  values  of  kl,  »«,  and  re  that  are  needed  in  the  hot-t^t  calculations.  For  example, 
Langdon^*  calculated  a  factor  of  two  reduction  in  kj,  doe  to  these  non-Maxwellian  distributions. 
Thus,  they  will  cause  both  upwards  and  sideways  shifts  in  the  Fig.  (3)  gain  curves. 

The  effects  of  generating  an  ion-acoustic  turbulence  at  the  critical  surface  can  also  significandy 
affect  the  gain  curves  of  Fig.  (3).  Some  of  the  more  important  effects  of  this  microturbulence 
might  be  modeled  through  the  introduction  of  a  single  parameter,  which  modifies  a  number  of 

important  plasma  properties.  For  example,  laser  absorption  would  be  increased:  =  ki,{1+/3). 
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and  the  electron-ion  coUisitm  time,  ^  =  '^6/  (1  +  0),  would  be  reduced  leading  to  inhibited 
electrcMi  heat  flow,  =  Q&/(1  -f  0)t  and  increased  coupling  between  the  electttm  and  ion 
thermal  energies  (see  Eq.  (4)).  Hig^  ion  temperatures  would  increase  die  amount  of  Doppler 
broactening  and  some^^iat  lower  the  gains  that  would  be  calculated  at  line  coiter  under  fairly 
mild  levels  of  saturated  ion-acoustic  turbulence.  Moreovn,  as  shown  in  Ref.  (18),  0  could 
acquire  values  between  1  and  5  at  die  plasma  densities  used  in  these  calculations  if  a  saturated 
microtuxbulence  amid  be  generated.  However,  determining  the  levels  of  microturbulence  that  can 
odually  be  generated  in  short-pulse  experiments  is  a  problem  currendy  uiufer  investigatitm.^’ 
Because  x-ray  laser  performance  is  a  sensitive  measure  of  underlying  plasma  conditions, 
die  observation  of  x-n^  laang  in  short-pulse  esqietiments  would  provide  important  tests  of  the 
dieoretical  models  that  ate  used  to  describe  laser  absorption,  transient  ionization,  and  x-iay  laser 
pumping.  The  peak  of  the  electrm  temperature  increases  in  much  the  same  wi^  as  the  absorbed 
laser  energy  as  a  functicm  of  (Hg.  (4)).  The  simultaneous  measurement  of  this  temperature 
alcmg  widi  the  gains  would  provide  useful  expmimental  information  on  the  presence  and 
influence  of  microturbulence  and  non-Maxwellian  electron  distributions  on  underlying  critical 
surface  absorption  jflieiuxnena. 
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Figure  Captions 

Figure  1.  Calculated  electron  (T*)  and  ion  (Tj)  temperatures  as  a  function  of  time  for  the  case 

=  2  X  10*®. 


Figure  2.  Calculated  average  charge  state,  {Z)  =  N^/Ni,  as  a  function  of  time  for  the  case, 
=  2  X  10*®. 

Figure  3.  Calculated  maximum  gains  for  the  two  J=0-1,  38-3p  transitimis  at  169  and  182  k 
versus  die  maximum  KrF  laser  intensity.  The  solid  J=0-1  (182  *A)  curve  is  corresponds  to  no 
modificatirm  of  die  J=0  m<»K^le  exdtadon  rate,  and  the  dashed  curve  corre^nds  to  a  faaor  of 
two  reduction  of  this  rate  (see  Ref.  (8)). 

Figure  4.  Calculated  absorbed  energy  per  micro-gram  of  sdenium  as  a  function  of  the  maximum 
KrF  laser  intensity  for  2  ps  FWHM  ganssian  pulses. 
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